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SPATIO-TEMPORAL EQUALISATION USING CHOLESKY FACTORISATION -AND SYSTOLIC ARRAYS 

5 

FIELD OF THE INVENTION 
The present invention relates to the field of the cellular telephone systems and in 
particular it relates to a space and time, quick convergence equalization method, for 

10 the cancellation of stationary and non stationary isofrequential interferers, in signals 
received from array antennas of base transceiver stations of a cellular mobile radio 
base system, and relevant systolic array equalizer. In particular, the invention refers to 
mobile systems employing a brief training sequence in transmission bursts. 

In telephone mobile systems, the main limitation to the traffic capacity offered to 

15 the public service is represented by the scarce band available. Therefore, all 
technological solutions enabling to increase the above mentioned capacity, or spectral 
efficiency, maintaining a given quality of the service in a wide area are particularly 
appreciated. To this purpose, appropriate cellular coverage strategies have been 
developed, together with multiplexing techniques in FDMA frequency (Frequency 

20 Division Multiple Access), and/or in time TDMA (Time Division Multiple Access), and 
more recently in space SDMA (Space Division Multiple Access) thanks to the usejDf 
the so-called intelligent array antenna combined to the receivers. 

The cellular coverage technique consists in dividing the service area in zones, 
called cells, each one served by its own base station BS (Base Station). Each cell 

25 avails of its own group of radio channels that can be reused in other duly spaced cells, 
in fact, the attenuation of radio signals is proportional to the distance with an exponent, 
or propagation coefficient, in general higher than 3. The whole of the cells employing 
the radio channels available to the system, without any reuse of the same, is called 
cluster of cells. Being the spectral efficiency inversely proportional to the cell area, 

30 microcellular systems are the most common approach to increase the spectral 
efficiency in high traffic density urban or suburban areas. As the cell area decreases, 
the problems associated to continuous handovers become more and more critical, 
consequently, mixed structures (microcells overlapped to umbrella macrocells) are 
generally preferred to simultaneously control quick mobile stations MS (Mobile Station) 

35 in scarcely populated areas, and slow mobile stations in high density areas. 
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In order to increase the spectral efficiency, a first approach as an alternative to 
the reduction of the cell dimension, consists in the reduction of the number of C cells 
per cluster, where C is also called cluster size, which is equivalent to increase the 
fraction of the total spectrum controlled by the single BS. This strategy is known as 
RCS (Reduced Cluster Size), whose theoretical limit is to render available the whole 
spectrum to all the cells, performing clusters of one cell only. 

A second approach consists in allocating a higher number of users in the 
spectrum already available to the cell, allocating more mobiles on the same channel, 
characterized by the same carrier and same time slot and therefore, without any 
intervention of the cluster size. This- second strategy is known as SCFR (Same Cell 
Frequency Reuse) and is based on the space division multiplexing SDMA within the 
cell. In practice, it is implemented in two different ways, the first one consisting in the 
division of the cell in smaller sectors through an arrangement of antennas with fixed 
radiation diagram (in general with 120° aperture), and a second and more advanced 
one, consisting in fitting the BS of relevant multibeam antennas, or synthetic aperture 
ones, that is capable to transmit and receive in the direction of mobile terminals with 
comparatively narrow directive beams; this implies that said beams must follow a 
generic mobile user moving in the azimuth plane. 

The formal difference between the use of an SCFR technique and the division of 
the cell itself in smaller sectors, resides in the capability to allocate the users in the 
relevant channels. Through sectorialization it is possible to manage more users on the 
same channel only if their position on the territory falls under different fixed sectors of 
radio coverage. In a scenario characterized by users in movement, this involves the 
need of continuous handovers. On the contrary, an SCFR structure is able to allocate 
more co-channel users in the same cell, without the need of partitioning the cell and 
have recourse to frequent handovers. Of course, even in this last case, the BS must be 
able to distinguish the co-channel users one from the other, in order to form 
appropriate radiation diagrams in receipt, each one capable to trace the relevant useful 
signal and reject the isofrequential interferers, corresponding to the remaining co- 
channel users simultaneously active within the same BS, and of course served by a 
same number of beamformers of the same BTS (Base Transceiver Station). As we 
shall see below, the discrimination between useful and co-channel interferers in 
GSM/DCS systems takes place availing of a training sequence made of 26 known 
symbols forming the midamble of the transmission bursts, assuming to have assigned 
different training sequences scarcely correlated among them to different co-channel 
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users. The above mentioned sequence is purposely foreseen for the extraction of the 
symbol synchronism for the estimate of the channel impulse response. In the case of 
co-channel users of the same cell, they are also synchronized among them, and 
therefore it is possible to locate the same through analysis of the training sequence 

5 only. 

When adopting the RCS and SCFR strategies, it is absolutely necessary to avoid 
that co-channel interferers can worsen the quality of the demodulated signal, thus 
making vain the advantage of the higher capacity of enabled traffic. We want also 
underline that, in the systems employing an approach of the RCS type, the interference 

10 is due to the co-channel users of adjacent cells, while in systems employing the SCRF 
approach the interference is generated both by the co-channel users of the same cell, 
and by the co-channel users of adjacent cells. 

The mobile radio channel has been studied for a long time, and different models 
have been proposed, consistently with the more frequent transmission environments, 

15 rural, urban, etc. From these studies it resulted that, due to the space and time 
dispersive nature of the radio channel, in addition to the interferers generated by co- 
channel users on the direct path, there are also those corresponding to the multiple 
paths. 

Particular attention shall be paid to possible delayed echoes of the useful signal, 
20 since in some examples of the known art they are confined to the role of co-channel 
interferers, while in the embodiment according to the present invention the same are 
actually recognized as integral part of the useful signal. 

Our attention shall be concentrated below, on base stations employing antenna 
arrays to implement the RCS or SCFR techniques, the latter also according to the more 
25 advanced SDMA technique, with dynamic assignment of the channels, or concurrently 
with the use of traditional techniques of fixed sharing of frequencies. 

We shall also start from the assumption that the synchronization problem has 
already been solved, that is the acquisition method of the carrier synchronism is known 
as well as the position of the training sequence inside the burst. 
30 BACKGROUND ART 

The nature of the mobile channel is such to cause a spatial dispersion (more 
arrival directions of echoes) and time dispersion (delay of echoes) of the power of the 
signal reaching each antenna. In a context characterized by propagation along multiple 
paths, it can be reasonably expected that the different paths of the signal reach an 
35 antenna array with different DOA. This characteristic of the mobile channel can be 
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used to perform an equalization in the spatial dimension, it is possible for instance to 
filter only the maximum power path, or appropriately combine the signals impinging 
from different paths among them to cancel co-channel interferes. The equalization of 
time dispersion shall be considered hereafter. 

5 In order to be able to effectively discriminate the useful from the different kinds of 

interfering signals, the array must be able to dynamically modify its own directivity 
function and insert some zeroes in its own radiation diagram in the directions of the 
interfering signals. We shall give below the description of some interfering implemented 
according to the known art, highlighting time by time the relevant drawbacks. The 

10 straightforward description, which however nothing detracts to the general character of 
the same, is limited to a linear and uniform array of antennas, called also ULA (Uniform 
Linear Array). 

Making reference to fig. 1 , the array consists of M equal aligned and equispaced 
elements. The separation between adjacent elements is d. For narrow band signals 
15 with frequency of the carrier fc, the corresponding wavelength can be defined as: 

J c 

where c is the propagation speed (c = 3.108 m/s in vacuum). The spacing d between 
the elements can be conveniently expressed according to the wavelength X as: 

rf = AxA (2) 
20 A is thus the separation between the elements in terms of wavelength. 

Assume later on, that a sole punctiform source is present, sufficiently distant to 
render the wave impinging on the array similar to a plane one. The direction indicated 
in fig. 1 by angle 9 is the direction of arrival of the plane wave, hereinafter called DOA 
(Direction Of Arrival). Therefore we obtain that the relative delay T between the signals 
25 received by two adjacent elements of the array is: 

dcos[9) 



30 



r = ■ 



(3) 



c 

the real signal received on the k-th element of the array is then: 

r k (t) = s{t-T k ) = s(t-kT) (4) 

It is assumed that the delays relevant to the different receivers along the array 
are small compared to the reciprocal of the signal band, so that the same can be 
approximated with opportune phase variations of the complex envelope r k (t) . 

Therefore it results: 
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r k {t) * s{t)e- J1 *< kT (5) 
the phase factor 2xf e r can be written in different ways: 

27tf,T = 2tf e = 2^co^6) = 2nA code) (6) 

C A 

Finally, if s{t) is the equivalent in base band of the signal received on the first 
element of the array, and defined a vector a(#) as: 

1 



e -J2x< A/-l)Acos(0) 



(7) 



we can group the signals received on the different elements of the array in a vector 

r(r): 

r(t) * z(0)s{t) (8) 
10 Vector a(0) is the so-called steering vector, or array response vector, at an 

impinging plane wave with DOA equal to 9 . 

In the background art, one of the first approaches to the problem of spatially 
selective reception and transmission is that of the so-called Switched Beams systems, 
that is those employing antenna arrays to form a fixed number of beams oriented in 
15 fixed directions. Beams so obtained, are selected, both in transmission and in 
reception, in order to employ the better oriented beam in the direction of the mobile 
terminal. Though easy to implement, these systems have the great disadvantage not to 
be able to effectively discriminate the interfering signals. 

Making reference to fig .2, we notice a Switched Beams receiver including a four- 
20 element array, connected to the inputs of a Beamformer block, whose four outputs are 
connected, to the inputs of a beam selection block and to those of a selector, 
respectively. The output of the beam selection block is connected to a control input of 
the selector. 

During the operation, the beamformer of fig. 2 is used to form a radiation diagram 

25 consisting of a fixed number B of beams in pre-set directions ^ on the azimuth plane. 
The inputs of the beamformer are represented by the vector r(r) (hereinafter denoted 
with y{t) ), neglecting for convenience, the tilde sign of the complex envelope. 

The beamformer consists of a matrix of coefficients [w mJ ,] defined as follows: 
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w 



m,b 



with m=1,2 M and b = 1 ( 2 B. The signal on the b-th output of the:beamformer 

shall therefore be equal to: 

y^S*) = f; j (/)e^' 2 ' T(m - |,4cos( ^ ) - cos(£,) ) (10) 

/7I= I 

5 if 9~(p b the output results the maximum one. In case more signals are present, 

impinging from several directions, considering the arrival direction 9-<f> b , the 
contributions of the different elements of the array are summed up in phase. The b-th 
output of the beamformer shall therefore be characterized by a radiation diagram with 
maximum in the direction <fi b . 

10 The selection of the best beam takes place according to a criterion tending to 

minimize the error probability, that is applying a MLSE estimate method (Maximum 
Likelihood Sequence Estimation) of the sequence transmitted at maximum likelihood . 
The above mentioned estimate is profitably made applying the well known Viterbi 
algorithm to a trellis sequential diagram, or trellis, whose branches represent all 

15 possible status transitions of the receiver. The algorithm calculates, for all the possible 

sequences {x n } , the distance between the corresponding signal transmitted: 

X'IM) = E„*.M/-*r) (n) 

and the received signal y{t) , and then selects the sequence [x n ] corresponding to 
the minimum one among said distances. From the above we can understand that the 
20 higher are the distances among the possible signals {*„}). the lesser is the 
probability to select the wrong sequence during the demodulation. Therefore, defining 
d min the minimum value of the distance among all possible y{t\ {x n }) , this parameter 

results being a reliable index of the error probability. 

Due to the calculation complexity of this approach, in some known examples 

25 another method is also proposed for the beam selection, simply based on the quality of 
the demodulation of the training sequence. Among all the possible beams, it is selected 
in fact the one determining the smallest number of errors in the midamble 
demodulation. At equal number of errors made, it is used the beam that, at the end of 
the Viterbi algorithm applied to the midamble, originates the minimum distance. 

30 An approach giving the antenna a higher intelligence level versus Switched 
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Beams methods, essentially consists in the execution of a spatial filtering on the 
signals received by the antennas, by a beamformer resuming in the output signal the 
directivity characteristics of an opportune radiation diagram capable to reduce the co- 
channel interference. This approach, in the major part of the known application 
5 examples, required the prior knowledge of the DOAs of the useful and co-channel 
interferes. Different types of algorithms studied to solve the problem of the DOA 
estimate are found in literature, availing of an array of receiving antennas. The most 
representative algorithm is that known as MUSIC (Multiple Signal Classification), 
proposed by R. O. Schmidt in the paper under the title "Multiple Emitter Location and 
10 Signal Parameter Estimation", published in IEEE Trans. Antennas and Propagation, 
vol. AP-34, No. 3, March 1986, reference to which shall be made for a more complete 
description. The algorithm requires the restrictive assumption that the covariance 
matrix of the signal is not singular; in this case the signals on antennas are not 
correlated among them and the steering vectors result linearly independent. 

15 When arrival directions are known / = 1, 9 p) an appropriate radiation 

diagram shall be constructed. A possibility is to define as many diagrams as are the 
estimated DOAs (p), on the basis of the MVDR technique (Minimum Variance 
Distortionless Response), described for instance by S. Haykin in the volume "Adaptive 
Filter Theory", published by Prentice-Hall, Inc., Englewood Cliffs, N.J., 1991. 
20 Said approach enables to estimate the beamformer coefficients in order that the 

radiation diagram has unitary response coinciding with the selected direction, 
minimizing at the same time the interference effect and of the noise coming from all the 
other directions. 

Therefore, up to this stage the arrival directions of the different signals have 
25 been estimated, making for each one of them an appropriate radiation diagram. 
However, according to the described technique it is not possible to know which is the 
arrival direction associated to the signal considered as useful and consequently, which 
is the radiation diagram more convenient to be used to perform the demodulation. One 
can think to select the beam according to the methods previously described, in order to 
30 select the one that enables a better demodulation of the training sequence. 

In the case of DOA estimate, the above arguments are strictly valid for the case 
of p signals uncorrelated with p<M f where M is the number of elements of the array, or 
for scenarios characterized by a low number of multiple paths, as the mobile 
communication systems in rural area. 
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In the case of propagation in environments characterized by number of multiple 
paths, such as for instance in metropolitan areas, the extraction of the DOAs can result 
ineffective, since the number of arrival directions of the different signals is often higher 
than the array size, and also the hypothesis of uncorrelation among the different 

5 signals on the antennas is scarcely verified. However, at least in the case where 
multipaths relevant to a same signal belong all to an angular sector, not too wide, the 
performances of the algorithm often remain satisfactory, and in any case higher than 
those guaranteed by the Switched Beams approach. 

To render the direction estimate more efficient in case of correlated signals, one 

10 can think to use other DOA research algorithms, like the WSF (Weighted Subspace 
Fitting) proposed by B. Ottersen and M. Viberg in the article "Analysis of subspace 
fitting based methods for sensor array processing", following the Proc. ICASSP 89, pp. 
2807-2810, Glasgow, Scotland, May 1989. This algorithm, at the cost of a higher 
computation complexity, is able to process also coherent sources. The advantages are 

15 however highly limited as already said both by the large number of arrival directions 
and by the fact that, even succeeding in estimating more than one direction for each 
signal, it should be necessary to use an algorithm for the beam composition capable to 
completely exploit the knowledge of said directions. . Availing of an effective 
demodulation diagram, capable to control the multipaths, as the one based on an 

20 MLSE algorithm, it can result useful to point the beam toward more replicas of the 
same signal, in order to maximize the signal-noise ratio, thing which is not possible with 
the already mentioned MVDR technique which, as already said, considers one sole 
direction as useful. 

The algorithms described above, foreseeing the DOA estimate prove also to be 
25 ineffective in case of non-stationarity during transmission. Let's consider in fact the 
case of an interferer disturbing the transmission only in a semiburst (and not during the 
' midamble). The arrival direction of an interferer, estimated using a beam selection 
algorithm based on the quality of the midamble demodulation, could not consider the 
effect of said interferer. In this case, it would be necessary to have recourse to more 
30 complex algorithms for the selection, or composition, of the beam based on the 
demodulation not only of the midamble but of the whole burst. 

In the cellular mobile system, the number of interferer signals may sometimes 
result higher than the number of the antennas forming the array, in this case the 
hypothesis p < M is no more verified and the array can also not be able to insert zeroes 
35 in the radiation diagram in the direction of each interferer. This means that the methods 
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15 



20 
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based only on the DOA estimate additionally loose their effectiveness. However, 
also in this context, the adaptive array could be able to combine the signals presents at 
the antennas to improve the ratio between the signal and the sum of interferes with the 
noise, ratio known in literature under the acronym SINR, (Signal to Noise and 
Interferers Ratio). In the following articles by J. H. Winters: 

• "Optimum Combining in Digital Mobile Radio with Cochannel Interference", 
published in IEEE J. Select. Areas Commun., vol. SAC-2, No. 4, July 1984; 

• "Optimum Combining for Indoor Radio Systems with Multiple Users", published in 
IEEE Trans. Commun., vol. COM-35, No. 11, November 1987, is developed the 
concept of optimum combiner. 

This is an alternative technique to the DOA estimate, which bases on the "prior" 

knowledge by the receiver of a sequence of known symbols {x n } (training sequence), 
present in the single transmission bursts for synchronization purposes. The sequence 
{x n } is additionally exploited in the spatial filtering process, through which the signals 
received from the antenna array are appropriately combined to maximize the SINR. 

This approach enables to avail of the spatial difference of the signal, besides 
reducing the relative power of the co-channel interferers, also to face the Rayleigh 
fading of the useful signal. 

Fig. 3 is the block diagram of an M elements optimum combiner. As it can be 
noticed, the combiner has M inputs for a same number of base band signals 
7i(«)».-.^/W»--^A/( w )- received by corresponding M elements of an array (not 
shown in the figure). The above mentioned signals are multiplied by the relevant 
weights w, ,...,vt>,.,..., w M , in general complex. Weighted signals reach the inputs of a 
first adder, at the output of which a y om («) signal is present, coming out from the 
combiner and reaches also, with negated sign, an input of a second adder, at the other 
input of which a x n signal is present. At the output of the second one, an error signal 
e n is present, reaching a WEIGHT ADJUSTMENT block having other inputs for signals 
y x (n) 9 ..- 9 y t {n) 9 .^ 9 y M (n). At the outputs of said block, weights w l ,.„,w / ,...,w M are 
present. 

In the operation, the signal y OM {n) can be expressed as: 



M 



(12) 
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with 



Denoting with w = [w l ,w 2 ,....,w /W ] r the vector of the weights and 

y (n) = [j/,(n), j> 2 1 y u( n )] 7 the vector of the signals received and' adopting a 
matrix notation, we can define the weighted sum of the signals coming from the 
antennas as: 

The difference between the desired signal x n , corresponding to the samples of 
the training sequence, and the signal y out {n) at the output of the combiner defines the 
error signal: 

e„=x„-w r y(/7) (14) 
10 This signal is used as control signal in the adaptation process of the combiner 

weights. 

The scope of the adaptation process is to find a vector of weights w, such to give 
at output, at any instant, a signal equal to the desired one x n . When the input signals 
can be seen as stationary stochastic variables, the WEIGHT ADJUSTMENT block 
15 searches the vector of the weights minimizing a cost function corresponding to the 
mean quadratic error MSE (Mean Square Error), defined as follows: 

r-(w) " E[\e n \ 2 ] (15) 

where £{A ] indicates the average value of the subject in brackets. 

By replacing the (14) in the (15) we can obtain the whole of the weights 
20 minimizing the MSE: 

e 2 (w) = E\\x\ 2 + w*y(«)y*(«)w-2w'x n y(«) 

1 _ _ 1 (16) 

where the spatial covariance matrix R^ of signals at the antenna input can be 
expressed as follows: 

R„ = £[y(«)y*W] (17) 
25 and the cross-covariance matrix R^ among the M signals at input of the antennas 
and the desired signal, can be expressed as: 

. = £[y("K] (18) 
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The mean square error defined in the (16) is a quadratic function of weights w. 
The components of the gradient of function MSE, are the partial derivatives of the MSE 
versus the value of weights. Obtaining the (16) versus w, we obtain the gradient as 
linear function of weights: 



= 2R w w-2R^ (19) 



When the selection of weights has been optimized, the gradient results null, then: 



R " W - = 1' (20) 

w ..„. = r:'r.„ 



opt 



The weight vector w opt is therefore the one giving the minimum mean quadratic 

error. The equation (20) is called Wiener-Hopf equation. For a more detailed analysis 
10 of the research methods of the optimum weight vector, reference shall be made to the 
article "Adaptive Antenna Systems by: B. Widrow, P. E. Mantey, L. J. Griffiths and 
BB. B. Goode, published in IEEE Proceedings, vol. 55, No. 12, December 1967. 

Considering now the time equalization, it is opportune to introduce the following 

considerations: let [x n ] a sequence of N symbols transmitted on a channel with 
15 dispersion length L, that is whose impulse response is significant within a time LT, 
being T the symbol time. Let's consider an array of antennas consisting of M elements. 

The vector of signals received on the array of antennas y{t) = [y^t), y 2 {t),....,y M {t)] 
is represented by the convolution of a vector h(r) without the transmitted 
sequence [x n } : 

20 yW= Z*» h ('-«r) + n(r) ( 21 ) 

where the vector elements h(r - nT) are the samples of the impulse response of the 
channel on the different antennas. The vector h(f) can be written as follows: 

h(0 = Z^^^o(^-r ( >te) (22) 

i = l 

being: 

25 • N the number of multiple paths, 

• A t the path loss of the i-th path, 

• 6. the arrival direction, 
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• h 0 the convolution of the modulation impulse and impulse response^ the reception 
filter, 

• a(6>.) is the steering vector, or vector of the array response in direction d t , defined 
in the (7). 

5 The vector n(r) represents the additional noise, we shall time and spatially 

consider white and that satisfies the following conditions: 

£[n(/V(0] = - 2 I (23) 
. *K (')«.(' + *)] = a 1 **) 
In the major part of the known realization examples of an optimum equalizer for 
the compensation of distortions introduced by the channel on time characteristics of the 

10 sequence transmitted, the equalizer is made up of a filter adapted to the channel 
impulse response, followed by a maximum likelihood detector (MLSE) based on the 
Viterbi algorithm. Both for the realization of the adapted filter, and for the MLSE 
detector one, it is absolutely necessary to known the channel impulse response, 
generally unknown. Availing of the knowledge of the midamble, the equalizer estimates 

15 the channel impulse response, to calculate then the coefficients of the filter matched to 
the channel impulse response and be able to estimate the sequence transmitted, 
profitably using the Viterbi algorithm. The possibility to perform an estimate of the 
vector h(0 makes therefore useless the explicit calculation of the same according to 
the (22), which presupposes the estimate of the arrival directions. 

20 Algorithms mentioned up to now, face the problem of the estimate of the 

sequence transmitted in two different equalization phases, namely: 

• spatial equalization for the cancellation of interferers; 

• time equalization of the spatially filtered signal. 

A separate equalization enables to spatially combine only the echoes of the 
25 signal resulting not significantly delayed. The delayed echoes of the signal are in fact 
considered as co-channel interferers. However, this approach can result in a reduction 
of the ratio between signal and noise at the time equalizer input, as it shall be noticed 
through comparison of figures 4 and 7, and in a lower robustness of the beamforming 
algorithm versus particular azimuth distributions of echoes, as it could be noticed 
30 comparing figures 5 and 8. 

In the light of what just said, let's consider the cost function (15) of the 
Wiener-Hopf combiner, previously dealt with: 
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e 2 (w I> w 2 ,....,vv (Vf ) = E 



= E 



M 



(15) 



where w. is the coefficient according to which the signal received at the i-th antenna is 
weighed. This algorithm maximizes the ratio between the signal and the noise, where 
signal means all those echoes resulting correlated with the training sequence {x n } . If 

5 the channel introduces delays higher than the symbol time, T~3J/js for GSM, 
echoes so delayed result decorrelated with the training sequence and are treated by 
the combiner as co-channel interferers. 

Fig. 4 is the radiation diagram on the azimuth plane of an antenna array, where 
the bearnforming coefficients are calculated with the Wiener-Hopf algorithm applied to 

10 the functional (15), in presence of an useful signal generated by a mobile MSd, with 
relevant echo delayed by 5 |ns, and of a co-channel interferer generated by a mobile 
MSL As it can be noticed, the radiation diagram shows a maximum in the useful 
direction (90°), a first zero in the co-channel interferer direction (120°), and a second 
zero in the useful signal echo direction (60°). Due to the second zero, the possibility to 

15 avail also of the power of the useful signal in the estimate process of the transmitted 
sequence is lost. 

Fig. 5 is the radiation diagram on the azimuth plane of the antenna array used in 
a scenario differing from the previous one only for the presence of a replica of the co- 
channel interferer placed in the same arrival direction (90°) of the main signal. As it 

20 can be noticed, the radiation diagram still shows the two zeros in the previous 
directions, of the interferer (120°) and of the echo of the useful signal (60°), 
respectively, but in this new scenario, the maximum in the useful direction (90°) is no 
more present. This behaviour can be explained, considering a traditional 
synchronization system, by the fact that the combiner sees a unique useful signal, 

25 different from the interferer, in the common direction (90°), however attempting to 
reduce the effect of the interferer in that direction reducing the gain. This involves an 
additional reduction of the power of the useful signal present in the estimate process of 
the transmitted sequence. Therefore the combiner looses the opportunity to transfer 
the maximum gain lobe in the sole echo direction of the useful signal (60°), thus 

30 succeeding to nullify the echo effect of the interferer having the same direction of the 
useful signal. 

In mobile radio environment more recent algorithms than those described up to 
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now are known, that try to overcome the drawbacks highlighted in the description of 
figures 4 and 5. The above mentioned algorithms still avail of the concept of 
minimization of the mean square error expressed by an appropriate functional, like it 
occurred for the optimum combiner, but, contrarily to this one, the information on the 
channel included in the expression of the above mentioned functional is now given by a 
so-called spatial-time filter due to the fact to adapt to the total response of the channel. 

An example of the above is represented by the international patent application 
WO 97 08849 A, published on 6/3/97 under the name of THOMSON CSF, under the 
title: "METHOD AND DEVICE FOR THE SPATIAL MULTIPLEXING- 
DEMULTIPLEXING OF RADIOELECTRICAL SIGNALS FOR SDMA MOBILE 
SYSTEMS", with French priority of 22-8-95. 

In the above mentioned application it is described a base station employing in an 
"intelligent" manner an array of antennas to which a multiway transceiver is coupled, 
equipped with beamforming means particularly suitable to enable the assignment of a 
same frequency-time channel to more mobiles in the same cell, provided that the 
relevant signals have sufficiently decorrelated arrival directions. The invention, which 
essentially claims an application of the SDMA technique in mobile radio environment, 
shows two embodiments distinguished by different typology of the beamforming means 
used in reception. More in particular, the first embodiment is referred to means 
requiring a functional block for the estimate of the steering vectors, which are then 
passed to a true beamforming spatial filter. The beamforming means of the first 
embodiment can be compared to those already mentioned for the DOA estimate and 
result far less adequate to the present invention compared to the beamforming means 
foreseen in the second embodiment, therefore the first embodiment shall not be 
examined. However, it is useful to remember that the techniques used in the first 
embodiment are called "transparent" because they do not require any prior knowledge 
of the sequence of symbols received. 

Beamforming means foreseen in the second embodiment of THOMSON paper 
use, on the contrary, the prior knowledge of the training sequence, as it occurs in the 
invention described below. These second means include a space-time filter. The 
optimization of the total response of the channel is made through two algorithms 
alternative between them, called FAST-BBT-DFE and FAST-DFE, respectively, where: 

• FAST is an acronym of Space-Time Matched Filter, 

• BBT is an acronym of Bruit Blanc Temporellement, with French terms, and 

• DFE of Decision Feedback Equalizer, with English terms. 
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We quote below what is said on the advantages of these algorithms: 
# -The FAST-BBT-DFE enables to rephase the power of all the paths of the useful 
signal and to reject in sub-optimal way the signals interfering from multiple paths. 
This procedure involves a gain in terms of signal to noise ratio on the current symbol 
5 compared to an only-spatial filtering algorithm (FAS-DFE), in fact the algorithm 
shows the same rejection capabilities of interferer signals with the difference that the 
echoes of the useful different from the main path are no more assimilated to 
interferer signals, so that a same number of freedom degrees are cleared to 
consider other interferers. 

10 • The FAST-DFE enables to process a signal with interferers coming from an array in 
an almost optimal way, since it performs not only the phasing of the echoes of the 
useful signal, but also that of interferers before their rejection. This procedure 
enables to eliminate an interferer source using one sole spatial freedom degree 
whichever is the number of echoes according to which this source is received by the 

15 array". 

However, it is necessary to point out that the above mentioned operations of 
"phase setting" or "re-alignment" are but the result of the concerned algorithms and do 
not represent autonomous steps of the channel optimization process, since otherwise it 
would require a separate estimate of the DOAs. Moreover, the DOA estimate block 
20 remains also in the second embodiment for the typical scopes of the SDMA technique, , 
that is to compile an intermobile correlation table of the arrival direction of signals to the 
antennas, to be used in the spatial multiplexing for the transmission beamforming. 

The algorithm FAST-BBT-DFE, described making reference to fig. 11 of the 
quotation, starts from the theoretical formulation of an optimum filter (being valid the 
25 restrictive hypothesis of the case) for each one of the N ways reaching the N antennas. 
The above mentioned filter does not directly act on the signal received, but on the one 
submitted before to a spatial whitening of the relevant noise+interferer component. To 
perform the optimum filtering it is necessary to known the noise+interferer 
autocorrelation matrix R b on the N antennas and the matrix G of total channel 
30 response. To this purpose, restrictive hypothesis are introduced on the nature of matrix 
that enable to separately estimate the two matrices and complete the filtering 
proposed at the beginning of the procedure. Then the filtered signal is allowed to reach 
a decider component of the DFE type. 

The main restrictive hypothesis consists in considering the matrix R b as having 
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impulse nature, which involves interferes completely uncorrelated between them, and 
therefore comparable to a white noise in the band of the useful signal. In other words, 
the interferes are considered free of echoes, or better, each echo is considered itself 
an autonomous interferer. Under these assumptions, the convergence of the algorithm 

5 appears rather quick. 

It is well clear that the application of such an algorithm is not recommended in 
mobile radio environment where, particularly in urban small size cells, the typical 
scenario is on the contrary characterized by co-channel interferes fit with a range of 
relevant echoes from multiple paths. The fact not to consider this reality involves the 

10 disadvantage not to be able to optimize in an actually credible way the total 
signal-noise ratio, or, whenever the SCFR technique (in SDMA mode) is used in order 
not to be aoie to separate the different co-channel mobile users in an optimum way on 
the azimuth plane. 

The FAST-DFE algorithm, shown making reference to fig. 12 of the quotation, 
15 does not place any limit on the noise nature, including the co-channel interferes, thus 
treating the channel in an almost optimal manner. Also in this case the structure 
consists of a space-time filter associated to a decision organ of the DFE type, even if 
the use of the Viterbi is considered possible, without describing methods for said use. 
The space-time filter consists of N time filters W, , N m long, one for each n-th way 
20 branching from the N antennas. An additional time filter HR is also present, it too N m 
long, placed on the feedback path of the DFE equalizer. As it is known from the theory 
of equalizers, a DFE, before submitting the symbols received to a threshold decider, 
tries to remove the ISI from the sequence of the relevant samples, providing to this 
purpose a deleter filter, acting in feedback on postcursor samples of the impulse 
25 response, and of a preliminary equalizer acting forward on precursor samples. In the 
diagram of fig. 12 the HR filter represents the deleter, while the whole of the filters 
whose filtered signals are summed up to form a unique sum signal, can be compared 
to the preliminary equalizer. The signal cleaned by ISI, denoted in the text with z{n) , is 
obtained from the difference between said sum signal and the signal filtered by HR. 
30 The z{n) signal crosses the threshold decider that gives an estimate %n) of the data 

sequence transmitted. 

The updating of the weights of the whole filtering structure is made by 
minimization of the mean square error MSE expressed by the following functional: 
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4 =E[|z(t)-3[t)| 2 ] (24) 
The functional (24) in case of stationary channel is estimated as follows: 

^—I 12(0-30! 2 (25) 

while the following expression applies to the non stationary channel: 

5 ? = - ± |z(i) - %if = X ?{n - 1) + |z(i) - i(0| 2 (26) 

where a is a forgetting factor selected between zero and one depending on the non 
stationariness of the channel. 

In presence of the training sequence, this replaces z(i) in the (25) and (26). 
The description of the second embodiment does not include a clear indication of 

10 how to act on the different filters to minimize the MSE t to description purposes, the 
acronyms of the known main minimization algorithms have been mentioned, but it is 
not supplied a valid example of the application to the specific space-time case. It does 
not also include the indication of a constraint to put on the power of the signal received, 
in order to avoid in the minimization of the (24) the convergence possibility toward a 

15 class of banal solutions. Assuming that in the THOMSON quotation, the embodiments 
with space- time filtering have been raised, to the rank of variant having its own 
originality, and that the application of the above mentioned algorithms does not seem 
known to the space-time case, the applicant is doubtful that the simple quotation of the 
same can enable those skilled in the field to reproduce the invention claimed in the 

20 embodiment without performing his own invention activity. 

ANALYSIS OF THE TECHNICAL PROBLEM 
What appears immediately clear is that a given complexity of the whole structure 
of fig. 12 of the THOMSON quotation. In fact, analyzing the expression of the functional 
to minimize, as well as the figure itself, and as clearly said in the text, we see that the 

25 optimization occurs both on the space-time filter W and on the HR filter of the DFE 
equalizer. Therefore, it is necessary to update a number of (JV + 1) x N HR independent 
coefficients at each symbol interval. The high number of coefficients to update involves 
a given slowness in the convergence of the estimate and therefore the need of long 
training sequences, this fact being autonomously highlighted in a table of the 

30 advantages and drawbacks of the FAST-DFE algorithm. 

If one wants to employ the above mentioned space-time equalizer in GSM 
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environment, where the training sequences foreseen by the standard strictly have a 
lengths of 26 symbols, and can therefore be considered short, the problem arises 
about the convergence of the estimate within the time slot of such a sequence. In this 
case, in fact, the estimate of error & would undergo an additional slackening due to 
5 lack of accuracy in the evaluation of the error itself, the DFE having to this purpose use 
the z{n) estimates in place of the no more available samples of the training sequence. 
Therefore, there is the actual possibility that a great part of symbols of the whole 
transmission burst is wrongly estimated, causing by this a degradation of the quality of 
the demodulated signal. Also, in case of transmission systems adopting short bursts, 
10 which could also include the GSM itself, with a burst of 148 symbols only, there is also 
the risk of a missed convergence within the whole duration of the same burst. 

The inconvenience highlighted is in itself sufficient to discourage the use of such 
an algorithm in GSM environment, but additional considerations shall be made to point 
out the consequent technical problems of non negligible importance. 
15 a second problem that can originate from the slowness in convergence of the 

estimate is that of the poor capability of the equalizer to face the channel non 
stationariness, mainly due to interferes that enable outside the time slot of the training 
sequence, and/or to the case in which the channel can no more be considered 
time-invariant, as it occurs for instance when the doppler effect is significant, due to the 
20 high speed of the mobiles. The above mentioned inconveniences become even more 
important in presence of sudden and high intensity interferes (c / 1 < -10 dB). This is 
a condition that can frequently occur in urban mobile radio environment, for instance, 
when due to a so-called "corner attenuation" the power of the interferes suddenly 
surpasses that of the useful signal. Particularly in this last context, more powerful 
25 solutions are required than that indicated in the (26), consisting in weighting each 
iteration of the estimate with the forgetting factor X < 1 . 

A last drawback, originated by the same causes of the previous ones, consists in 
the difficulty to physically realize a receiver able to perform in one sole symbol time all 
the arithmetic operations necessary to the adaptation of a so great number of 
30 coefficients of the space-time filtering structure. 

SUMMARY OF THE INVENTION 
Scope of the present invention is to overcome the above mentioned drawbacks 
having as main scope that to indicate a spatial -time equalization procedure particularly 
quick in converging towards optimum performances, and that can be profitably 
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employed in base stations of a cellular mobile system, for instance GSM, fit with so- 
called smart antennas. 

A second scope, concurrent with the first one, is to indicate the method to extend 
the equalization procedure outside the midamble, in presence of a non-stationary 
5 channel, particularly if it is affected by sudden and high intensity interferers, assuring 
the shortest possible duration of the transient within which a loss in quality of the 
demodulated signal occurs and the limitation of such a loss to limited values. 

A third scope is to indicate an iterative digital method for the updating of an 
estimate of the optimum equalization parameters, accurate and of quick execution at 
10 the same time. 

Fourth scope of the invention is to indicate the realization methods through 
systolic arrays of an equalizer based on the just mentioned iterative digital method. 

To attain these objects, scope of the present invention is an equalization method 
to cancel the isofrequential interferers in signals received by array antennas of base 
15 stations of a mobile radio system. Said signals being received as bursts of symbols of a 
transmission burst, including a training sequence known to the receiver, said received 
bursts being demodulated, digitized, stored, and synchronized through discrimination 
of the training sequence. 

The main technical problem highlighted in the known art is solved by the 
20 equalization method of the present invention, through the following operational phases: 

a) spatial filtering of signals coming from antennas, using a spatial filter combining 
corresponding symbols among them on the basis of appropriate weights in the 
different bursts stored; 

b) time filtering of the signal generated by said spatial filter, using a filter adapted to 
25 samples of the impulse response of the transmission channel; 

c) combined optimization of said weights and samples on the basis of the symbols of 
said training sequence and said signals coming from the antennas present in bursts 
stored; 

d) estimate of symbols of the transmission burst corresponding to the sequence 
30 coming out from the time adapted filter; 

e) updating of the optimized values of said weights and samples, out of the time slot of 
the training sequence, through repetition of the previous phases a)...d) ( replacing in 
phase c) the symbols of the training sequence by the symbols estimated in phase 

• d), thus rendering said equalization adaptive, as described also in the first claim. 
35 Profitably the decider component is a Viterbi processor exploiting said optimized 
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samples of the impulse response of the channel to calculate, in the known way, branch 
metrics of a trellis, used to estimate at maximum likelihood the sequence of symbols 
transmitted. The decider could also be a DFE. but this would involve a loss in 
performances. 

Said combined optimization of the equalization parameters is obtained through 
minimization of a functional corresponding to the mean quadratic error originating from 
the comparison between corresponding w-th samples of two sequences, of which, a 
first one is obtained submitting the symbols of the training sequence to a time filter that 
estimates said samples of the impulse response of the transmission channel, and a 
second sequence is obtained submitting the «-t samples of said stored bursts to a 
spatial filter estimating said weights. 

Profitably, the above mentioned combined estimate of the optimum equalization 
parameters is made in an iterative way on complex numerical magnitudes expressed in 
algebraic matrix form and representing the envelopes of relevant signals in base band. 
15 Compared to the more pertinent background art consisting of the FAST-DFE 

method, the equalization method scope of the present invention has the advantage of 
a quicker convergence, originating the capability to better cope with particularly critical 
transmission scenarios. This advantage proves to be essential in GSM environment 
and in systems adopting particularly short transmission bursts. 
20 The advantages highlighted essentially derive from the fact to perform the total 

filtering in two subsequent phases, a first only spatial phase, and a second only time 
phase, but to optimize at the same time the coefficients of both the filters (as for the 
coefficients of the time matched filter, the same are indirectly optimized through 
optimization of the samples of the impulse response). Thus a reciprocal iteration 
25 determines between the two types of filtering tending to the optimum due to the effect 
of the constraint imposed by the MMSE criterion (Minimun-MSE) simultaneously to the 
sure data of the training sequence. The result of such an interaction is to aid the 
crossing through the cascade of the two filters of the maximum possible power coming 
from the useful signal, therefore including its echoes and to reject the major possible 
30 number of co-channel interferes, to which the relevant echoes must be related. What 
is actually made, is an automatic resynchronization of echoes, valid both for the useful 
and for the interferes, which operates in synergy with the discrimination of the direction 
of the different signals- on the azimuth plane. The functional characteristics of the 
subject procedure distinguish it in an original manner compared to the procedures of 
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the more pertinent known art, where, as we have seen, foresee one sole space- time 
filtering phase to completely shape the channel, which can be followed by an 
additional time-oniy filtering phase in case a DFE is used. 

Thanks to its functional characteristics, the procedure according to the invention 
5 handles a considerably limited number of coefficients, equal to the sum M+N of the 
coefficients of two FIR filters, this undoubtedly facilitates the estimate convergence 
speed. 

Additional advantages relevant to the calculation methods and to the non- 
stationary channel shall be listed time by time as the description goes on. 

10 It must be pointed out that, notwithstanding the similarity of the form of the 

functional (26) and that used in the present invention, which shall be described below 
and a clear representation of which has already been given, actually there is a 
difference in the two terms of the difference forming the subject of the module. In the 
case of the known art, in fact, the correction action of the second time-only term is less 

15 effective in modifying in the optimum direction the spatial component, compared to the 
action of the corresponding term in the functional of the present invention, because a 
spatial-only component does not exist in the reality, being the first term of space-time 
nature, and this is an additional cause slackening the convergence. 

Finally, the term "joint estimate" in the present application assumes an effective 

20 meaning supported by the description, distinguishing by this from an implied meaning 
that the same wording has in the description of the FAST-DFE algorithm. The 
execution method of such an estimate shall be given in relation to the description of an 
innovative method for the minimization of the functional error, which shall be the scope 
of relevant ensuing claims. It is worth to mention that said method is ground on the 

25 iterative updating of a matrix structured by an updating vector including, the signals of 
the stored bursts and the symbols of the training sequence, and of the use of the 
updated matrix for the calculation of said weights and samples optimized at the n-th 
symbol time. 

In the method main scope of the invention the combined estimate of the best 
30 equalization parameters exhausts at the last sample of the training sequence, generally 
included in a midamble of the transmission burst. The 26 samples of the GSM 
midamble (out of the 148 available) are sufficient to assure the convergence of the 
estimate even in presence of medium quality signals received. The best equalization 
parameters remain then available for the demodulation of the remaining part of the 
35 stored burst, and the demodulated signal does not loose quality, provided that 
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channels are time-invariant and the interferes have stationary characteristics. 

The method according to the invention solves in the adaptive phase e) the 

T 

problem of a correct equalization though in presence of channels with time-variant 
characteristics and/or non-stationary interferes, updating the above mentioned 
5 estimate of the best equalization parameters at each new symbol received. From the 
operational point of view, it is necessary to this purpose to avail both of a vector of 
antenna input signals, and of a vector of data transmitted. While the first vector is 
always available, the second vector is completely known only during the midamble, 
where it corresponds to the training sequence stored in the receiver. We see therefore 
10 the need, once the midamble is exhausted, to replace the training sequence by the 
symbols of the estimated sequence, considering that the symbols already decided at 
an instant n, are but those associated to a path selected in the trellis up to instant n-th. 
Simultaneously, an updated sequence of optimum samples of the impulse response of 
the channel is sent to the Viterbi processor, for the updating of the values of the branch 
15 metrics. With this operation, the advantages mentioned above can be completely 
appreciated along the whole transmission burst. 

This way of proceeding inevitably introduces an excessive delay in the estimate 
of the optimum equalization parameters, due to the reaction time of the Viterbi 
algorithm which adds to the delay due to the time matched filter, thus preventing to 
20 follow quick variations of the channel characteristics. 

The inconvenience is avoided forcing the Viterbi algorithm to decide just after few 
steps (typically 4, or 5) which is the most possible sequence of data, that is taken at the 
time of the present symbol climbing in the trellis the maximum likelihood path, while 
final decisions for the MLSE detector are taken at the end of the burst, in order that the 
25 same are as much reliable as possible. 

Since short term decisions can result wrong, an error propagation phenomenon 
can be generated, which can last until wrong data are reused for the estimate of the 
best equalization parameters, that is for a number of symbol times equal to the intrinsic 
delay the short term decisions are made available. 
30 The equalization method scope of the invention, eliminates also this second 

drawback copying as a whole, at any symbol time, the L symbols necessary to the 
updating of the equalization parameters. It is well evident that the updating of the best 
equalization parameters is but a new combined estimate of the same. 

Profitably, it is possible to proceed to the above mentioned updating, and to the 
35 consequent adaptive equalization, starting from the ends of the midamble and 
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proceeding in parallel in the two directions towards the ends of the burst previously 
stored and synchronized. This involves the advantage to avail of the double of the time 
for calculations. It is also possible to proceed with a first pass on the right semiburst 
and on the left one, only afterwards, using the decided data of the right semiburst to 

5 have a longer group of data on the left one. This procedure improves the accuracy but 
not the speed of the estimate. 

With the actions made up to now, and considering that the combined estimate of 
the optimum equalization parameters is extended to the whole burst, it is useful to 
increase the effectiveness of the more recent data versus older ones, in order to render 

10 the method more reactive in respect with the variations of the channel characteristics. 
To the purpose, it is introduced a forgetting factor lower than one unit, as indicated on 
the other hand also in the known art, which appropriately weighs the elements of some 
matrices storing the historical data used in the estimate, before any updating iteration 
of the same. 

15 The equalization degree obtained with the method described is generally 

sufficient to enable a correct demodulation in almost all the real transmission 
situations. However, in presence of high-intensity non-stationary interferes 

(c//<-10dB), wrong decisions taken during the starting phase of the adaptation, 

immediately after the arrival of the noise, are so numerous to overcome the threshold 

20 enabling the estimate of the best equalization parameters to converge towards BER 
acceptable values in the demodulated signal, in a reasonably short time versus the 
semiburst length. The high increase of the error probability is unavoidable due to the 
intrinsic delay of short term decisions, which slackens the instant of interferer view, and 
that of the implementation of the countermeasures to neutralize the same. 

25 To face high intensity non-stationary interferes, it is described an embodiment 

including a phase of continuous monitoring of the power received on the antennas, and 
of non-coherent detection of possible peaks higher than a given threshold. This phase 
is profitably performed in parallel to the normal operation. The detection of a peak at an 
instant n is used to automatically select the operation according to the embodiment. 

30 After that, for the whole AT duration(typically twenty samples) of a convergence 
transient of the above mentioned matrices of historical data used in the estimate, only 
the coefficients of the spatial filter are updated and not the estimated samples of the 
channel impulse response, which maintain the values they had at instant n . This 
strategy is aimed at anticipating at instant n the future knowledge on the interferer to 
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cancel, present in the samples of the stored burst. In this way the interferer is "seen" in 
advance in the calculations and the adaptive error minimization algorithm can start 
updating the beamformer coefficients to try to eliminate the noise. Of course such an 
advance is not possible for the samples of the impulse response, whose adaptive 
5 estimate mandatorily requires short term decided data. 

The estimate as in the embodiment results therefore "vitiated", because to obtain 
the advantages it aims at, it introduces a synchronization loss between the samples of 
the signals coming from the antennas and the relevant decided data. The choice to 
maintain the estimated samples of the channel impulse response constant, reveals 
10 therefore as a caution measure destined to avoid error propagation, otherwise due to 
the sudden updating of the sole spatial component of the structured matrices of 
statistical nature. 

The original operation is automatically re-selected when the transient of AT 
duration elapses. 

15 The advantage obtained by the embodiment is to render the convergence of the 

equalization parameters estimate more effective, also in presence of sudden and 
strong interferers enabling out of the midamble, reaching at the end of the brief 
transient a more robust cancellation of the echoes of the interferer. 

Therefore another scope of the invention is a method embodiment, implementing 

20 the teachings just described, as described also in the relevant claims. 

It is worth highlighting that the processing of the interferer as indicated in the 
embodiment, which requires to have separate matrices for the space and time 
component of the signal, cannot be conceived with space-time filtering structures like 
those of the known art. 

25 To make the point on the sequence of equalization phases implemented by the 

subject method for the demodulation of the whole burst, the following is obtained: 
equalization starts in presence of the training sequence, once the midamble is 
exhausted it proceeds using low delay estimated data, the possibility to automatically 
enable in parallel the operation according to the alternative embodiment should be 

30 provided. Should this occur, the alternative embodiment quickly exhausts its effects 
within the time AT , at the end of which it is convenient to automatically reselect the 
previous operation, all within the time interval of the 0,577 ms burst. 

The equalizer based on the subject method, and its alternative embodiment, 
does not simply repropose a device described in a functional way retracing the steps of 

35 the method, but it exploits the capabilities of systolic array implemented through digital 
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integrated circuits of the ASIC type (Application Specific integrated Circuit) purposely 
configured. 

Consequently, additional objects of the present invention are equalizers 
operating according to the equalization method scope of the invention, and its 
5 alternative embodiment, and include systolic arrays for the combined estimate of the 
optimum equalization parameters, as better described in the relevant claims. 

BRIEF DESCRIPTION OF DRAWINGS 
The many advantages and scopes of the present invention will be appreciated 
from the following detailed description of a preferred embodiment of the same with 
10 reference :o the attached drawings, in which: 

- Fig 1 snows an array of antennas (array) used in the present invention; 

- Figures 2 and 3 show the block diagrams of a same number of beamforrners made 
according to the prior art; 

- Figures 4 and 5 are the radiation diagrams on the azimuth plane, obtained applying 
15 the beamformer of fig. 3 to an array similar to that of fig. 1 , concerning two different 

transmission scenarios in mobile radio environment; 

- Fig. 6 is a general block diagram of a receiver particularly emphasizing an equalizer 
operating according to the equalization method of the present invention; 

- Figures 7 and 8 are two radiation diagrams on the azimuth plane, obtained applying 
20 to the array of fig. 1 a beamformer included in the equalizer of fig. 6, relating to two 

transmission scenarios identical to those that originated the radiation diagrams of 
figures 4 and 5; 

- Fig. 9 shows two time successions of images placed side by side, ordered from top 
to bottom for growing values of the discrete instant n, out of which, the left 

25 succession indicates the radiation diagrams on the azimuth plane obtained applying 
the beamformer of fig. 6 to the array of fig. 1 , and the right one shows the estimated 
channel impulse response; 

- figures 10 and 11 show the BER curves at the output of the receiver of fig. 6, 
according to the possible length of the training sequence of the transmission burst, 

30 for two different values of C//; 

- Fig. 12 shows a general block diagram of a receiver particularly highlighting an 
equalizer operating according to the equalization method of the present invention 
and its alternative embodiment for strong interferers; 

- Fig. 13 shows a trellis relevant to a STIM block included in equalizers of figures 6 
35 and 12; 
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- Fig. 14 shows two trends of BER along the burst, evaluated at the outputs of 
receivers in figure 6 and 12, respectively; 

- Figures 15, 16, and 17 show the radiation diagrams on the azimuth plane relevant to 
the curve shown with solid line in fig. 14; 

5 - Fig. 18 shows three trends of BER along the burst, evaluated at the output of 
receiver of fig. 12, for three different values of the C/l ratio of the input signal; 

- Fig. 19 shows to comparison purposes two trends of BER along the burst, evaluated 
at the output of receiver of fig. 12, concerning the main method and its alternative 
embodiment; 

10 - Fig. 20 is a block diagram of an estimate block FUNZ included in equalizers of 
figures 6 and 12; 

- Fig. 21 is a block representation of a first systolic array included in the FUNZ block 
of fig. 20; 

- fig. 22 is a block representation of a second systolic array included in the FUNZ 
15 block of fig. 20; 

- fig. 23 shows a class of additional blocks included in the second systolic array of fig. 
22; and 

- fig. 24 shows the signals involving the generic blocks of the systolic array of figures 
21 and 22. 

20 DETAILED DESCRIPTION 

Making reference to fig. 6, we notice the general block diagram of a receiver 
adopting an equalizer operating according to the method of the invention. The receiver 
is connected to an array ARY of M antennas giving signals, represented as a whole 
with an algebraic vector y{n) . The vector y{n) reaches respectively the inputs of a 

25 block BMF and the first inputs of a block FUNZ. The BMF block is a spatial filter, or 
beamformer, reached by a vector of the optimized weights w op ,(«), coming from first 
outputs of the FUNZ block. At the output of the BMF block there are the real and 
imaginary (l,Q) components of a signal y out {n) which is filtered by a FADAT filter 
matched to the impulse response of the transmission channel. The components of the 

30 base band signal filtered by FADAT reach the first inputs of an estimator STIM, at the 

output of which an estimated sequence {B n } is present. The second input of the FUNZ 

block is reached by a training sequence {x n } stored in the receiver, previously 
synchronized to a copy of the same present in the midamble of the transmission burst, 
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availing to the purpose, known methods based on the oversampiing of the signal 
received and the correlation of the sequences generated by the oversampiing with the 
sequence stored in the receiver. A vector of L optimized samples h opt {n) of the 

channel impulse response comes out from second outputs of the FUNZ block, directed 
5 towards additional inputs of the FADAT filter and towards second inputs of the STIM 
block. 

The FUNZ block includes in its turn a PESI block that calculates a vector of M 
weights w(«) of the beamformer BMF, and a CAMP block that calculates a vector 
h(n) of L elements, corresponding to a same number of significant samples of the 

10 impulse response of the transmission channel. It includes also a two-input adder 
^ reached by two digital signals obtained from appropriate processing of PES! and 
CAMP blocks, from which the signal coming from CAMP is taken with sign minus. At 
the output of the adder an error signal e n is present that returns to the above 
mentioned blocks. The first and second inputs of the FUNZ block correspond therefore 

15 to the inputs of PESI and CAMP filters, respectively, the outputs of the same filters are 
also the first and second outputs of the above mentioned block. 

The representation of the FUNZ block highlights only the main aspects of the 
calculation of w opf (w) and h opt {n) t a more detailed description of the internal 

architecture of this block shall be made describing fig. 20. 

20 For ease of discussion, the diagram of fig. 6 does not include some blocks 

referred to known operations made by the receiver and necessary to the operation of 
the equalizer, such as for instance those relevant to the carrier synchronization, the 
conversion/demodulation in base band of the M bursts received from a same number 
of antenna, the sampling and analog/digital conversion of the same, the storing of the 

25 M bursts of digitized symbols, the acknowledgement of the training sequence and the 
synchronization of stored bursts, etc. Furthermore, all the elements of vectors are 
intended referred to complex envelopes in base band. In particular, the y(n) signal 

coming from the antenna is that defined in the (21). 

In the operation, the spatial filter BMF is a MTAPS transversal filter obtaining the 
30 signal y out {n) at its output making the weighed sum of the y(n) signals coming from 
the antennas and present in the bursts stored: 

^W = w r yW (13) 
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The adapted filter FADAT and the estimator STIM require, for their correct 
operation, to know the vector h 0/M (w), after that their realization is known to the 

technician of the sector. In particular, the FADAT filter is a LTAPS transversal filter, in 
genera! with frequency response of root raised cosine type, where L is the channel 
length. On the contrary, the STIM block is a processor performing an MLSE estimate of 
the transmitted sequence on the basis of the well known Viterbi algorithm, applied to a 
trellis with branch metrics whose expression is the known one of the family of receivers 
operating with pilot symbols (training sequence) and differential coding. 

In the equalization method governing the operation of the equalizer of fig. 6, the 
estimate of vectors w op# («) and h opt {n) is made jointly, considering the following cost 

function. 

* ,2 {"i , »\vm A-i) = E k 



= £ 



; = 0 



(27) 



from which it results that the mean square error e 2 , differently from the Wiener-Hopf 
optimum combiner, results a function both of coefficients w. of the beamformer and of 
the samples of the channel impulse response h t . . In this case the time dimension is 
directly introduced in the error functional to minimize. 

It is now proposed a matrix notation, to the purpose to render the processing of 
the problem of the minimum of the (27) more compact. Denote with Y the matrix 
corresponding to N time samples at the M antennas: 

y r (») 
y r U-D 

M 

_y r («-^ + l)_ 

obtained from vector y(n) = [y x (n), y 2 («), y M (n)] T whose M elements are the 

samples of signals received on the antennas; therefore: 



Y = 



(28) 



Y = 





y n +\ 


A 


y n *M-\ 


y n -\ 


y n 


A 


y n *M-i 


M 


0 


O 


M 




y n -N+2 


A 


y n -N+M . 



(29) 



25 



Concerning the number N, we shall see in the description of figures 10 and 11 
that a value included between 20 and 30 ( N = 26 for the GSM midamble) is sufficient 
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to furnish a good statistical basis for the estimate of the best equalization parameters. 

Denote with X the matrix of data transmitted (limited just for now to the training 
sequence), each one having duration L: 

x r (") 
M 

x T (n- N + 1) 



X = 



(30) 



obtained from the vector x(n) = [.v„("),-v„ +1 («),..:.,.v„, i _,(«)] r whose L elements are 
data transmitted coinciding with the training sequence \x n ) , and duly stored, assuming 
a channel length equal to L; so that: 



X = 



X„ X n^ 


A 




■Vl X n 


A 


X n-rL-2 


M O 


o 


M 


X n-N + \ X n-M-r2 


A 


X n-H+L_ 



(31) 



and with: 

10 w = [w 1 ,>v 2) ....,w A/ ] r (32) 

h = [h l9 h 2 ,.... 9 h L ] T (33) 
the vectors of the coefficients of the BMF spatial filter and of the channel impulse 
response are respectively indicated. 

The problem proposed ip the equation (27) can therefore be rewritten in the 
15 matrix form, as follows: 

(w op ,,h op ,) = argmine 2 (w,h) = arg min||Yw - Xh| 2 (34) 

||h|| 2 = 1 (35) 
where in the (35) is introduced the constraint on the power of the impulse response to 
avoid the banal solution (w op , =0,h op , =o) . The constraint (35) is a consequence of 

20 the fact that the bits of the sequence transmitted are assumed as time uncorrelate, 
being valid other hypothesis the (35) should be modified in a consistent way. 

The problem of the minimum constrained can be faced applying the Lagrange 
method. To the purpose, denoted with y a real parameter (Lagrange multiplier), let's 
consider the following functional: 
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j(w,h,A) = || Yw - Xh|| 2 + r(|hf " l) 06) 
The research of the minimum constrained in the (34) reduces -then to the 
research of a "solution" (w, h) of the problem of free mrnimum of the functional 
y(w,h,;/) . The above mentioned solution can be searched considering the system: 

dJ_ 



= y w (w,h,v) = 2Y"(Yw-Xh) = 0 

= y h (w,h,/) = -2X*(Yw-Xh) + 2>h = 0 



(37) 
(38) 



J£ = y,(w,h,y) = h'h-1 



(39) 



The equation (37) can be construed as the introduction of the error orthogonality 
principle, in fact it is possible to express the equation (37) as: 

E[e n y'(n)]=Q (40) 

Resuming, for convenience, some expressions partly already introduced, denote 
with R w the spatial covariance matrix expressed by: 



R„ = E[y{t)y(t)l 



(41) 



and likewise with and R^ , respectively, the covariance matrices on the data and 
of cross-covariance between data and inputs to the antennas expressed as: 



def - , 

R„ = E[x{t)x{t)\ 

R 9 = E{x{t)y(t)\ 
We shall use below the sampling co-variance matrices defined as follows: 

1 N 



1 N 



n-l 



We can therefore define a matrix R structured as follows: 

R = 



Ryy 



R 



(42) 
(43) 

(17) 
(44) 
(45) 

(46) 
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which shall enable a more simple representation of the problem. Equations (37) and 
(38) can be re-written respectively as: 

R >yW = R'^h (47) 

R^w-R^h = yh (48) 

From the equation (47) it is possible to obtain w as function of h and replace it 

in the (27) to obtain the mean square error. While replacing the w value in the (48) the 
following equality can be demonstrated: 

? : (w,h) = h'R7h = -v (49) 

where: 



is called Schur complement of R vy in R. The estimate of vector h(n) that minimizes 
the mean square error, corresponds therefore to the estimate of the eigenvector 
h in) of matrix R~ associated to the minimum module eigenvalue y . 

Opt x * 5 

Therefore it is possible to estimate also the w c/w («) coefficients of the 
15 beamformer replacing vector h opt (n) in (47): 

Summarizing, the choice of the optimum weights minimizing the cost function, 
can be synthesised in the following main points: 



1. calculation of R^, R^ matrices and of the Schur complement R, ; 



2. estimate of the eigenvector of matrix R^" associated to the minimum module 
eigenvalue h opt (n)\ 

3. estimate of w («) coefficients of the beamformer. 

Assuming to implement the equalization method summarized under points 1., 2., 
3., with a program controlling a finite precision processor, the calculation of the Schur 
25 complement Rf can result a very critical operation, mainly if the word length of the 

processor is very limited. 

The above mentioned criticality is considerably reduced thanks to the use of the 
methods developed in a theory on partial triangular factorization of the matrices, known 
in literature as "Displacement Theory", or even "Displacement Structure", included in 
30 the bibliography given below, reference to which shall be made. The key element of 
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this theory is an entity called "Cholesky factor" of a matrix, in the case considered 
matrix R, whose calculation is made through the simple updating of a previous value, 
intervening with appropriate transformations (Givens rotations) directly on the new 
data. As we shall see, the updating of the "Cholesky factor" gives a matrix, structurally 

5 different from the R but having the same information content, whose elements highly 
facilitate the w opl and h opl calculation. It is thus excluded the explicit calculation of 
sampling covariance matrices and the recourse to critical operations, such as that of 
the direct matrix inversion, since the above mentioned rotations do not involve any 
calculation of the determinant. This procedure involves a better conditioning of the 

0 calculation problem, and a consequent reduction of the minimum word length of the 



Such an approach gives the additional advantage to extend the calculation to the 
whole burst; the solution described is in fact intrinsically recursive (iterative) . On the 
other hand, a non recursive method, that is involving the explicit calculation of the 
15 sampling covariance matrices, is characterized by a high complexity in the calculation 
of the Schur complement through the (50), which makes the extension of such a 

method to the whole burst difficult. 

It must be underlined that the application examples of the "Displacement Theory" 
known up to now, concern only the optimization of the performances of time-only 

20 equalizers. A new element has been that to connect the mathematical formulation of 
the solution of the problem of minimum of the error functional used in the present 
invention, to the formulation given by such a theory for the solution of the problem of 
minimum in MMSE-DFE time-only equalizers. Due to the similarity highlighted between 
the two formulations it has been therefore possible to avail of all the knowledge 

25 developed for time equalizers adjusting it to the space and time equalizer of the 



This stated, let's consider the matrix R described in the (46) which, being 
positive defined, assures the existence of a unique Cholesky factor L such that: 



One of the fundamental methods for the triangular factorization, is the Schur 
reduction process described: 

• in the article by J. Chun and T. Kailath, under the title "Generalized Displacement 
Structure for Block-Toeplitz, Toeplitz-Block, and Toeplitz-derived Matrices", 



processor. 



invention. 



R = LF 



(52) 
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published in NATO ASI Series, vol. F 70; 
• in the volume "Numerical Linear Algebra Digital Signal Processing and Parallel 
Algorithms", authors G. H. Golub and P Van Dooren, edited by Springer-Verlag, 
1991. 

The process enables, after M factorization steps, the reduction of matrix R as 
follows: 

"0 0 
0 R. x 



| L 

= L U J 



[L- IT] + 



(53) 



where can be then considered as the residue of the partial displacement in 
triangular form. This approach enables to calculate the matrix without having to 
10 invert the spatial covariance matrix R xy . Starting from the Schur complement, obtained 
from the factorization process, it is possible to estimate h opt according to classical 
estimation methods of the eigenvector associated to the minimum module eigenvalue, 
for instance, the method of inverse powers. 

The partial factorization highlights also the possibility to have matrices and 

15 R ^ factorized in their triangular form: 

R„ = LL* (54) 
R^ = UL* 1 (55) 

enabling the reduction of the equation (51 ) in the form: 

L'w 0 „ = U'lv (56) 
20 The lower triangular form of matrix L f as we shall see below, enables the 

solution of the equation (56) with a simple process called "back-substitution" process. 

Let's now analyse the different sequential steps of the equalization method 
employing the Cholesky factor in the calculations, giving for each of then the relevant 
complexity: 

25 • step 1: application of the known Schur algorithm for the quick reduction of the 
matrix R for the Cholesky factorization; 
• step 2: application of an algorithm known in the literature as "Method of inverse 

powers" for the estimate of the eigenvector (h op/ ) of associated to the minimum 

module eigenvalue; 
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• step 3: solution of the equation (56) for the estimate of optimum coefficients of the 
beamformer (w opf ), according to a method known in literature as "back-substitution" 

method. 
Step 1 

5 From the Toeplitz theory of matrices it is known that a generic matrix R{n) has a 

displacement structure Toeplitz-type if the displacement matrix VR(«) defined as: 

VR(n) = R(n) - F(«)R(w - l)F'(n) (57) 
has low rank (rank(VR)"rank(R) ). The matrix R expressed in the equation (46) falls 
within this category and, in accordance with the above mentioned theory, it can be 
10 calculated in an iterative way, arbitrarily assuming R(0) = 0, and extending to the 
whole midamble the updating defined as follows: 

R(n) = R(»-l) + g„g; (58) 
where the updating vector g n , denoted as: 

"y(») 



gn = 



x(«) 



(59) 



15 is called generator of the difference matrix, since imposing F(n) = I we can write 
VR(«) = g n z' n . 

Likewise, it is possible to update the Cholesky factor of matrix R{n), thanks to 
the special time-varying nature of the same. Therefore, denoted with L(«) the 
Cholesky factor of matrix R(«), such that R(«) = L(h)L*(w) , it is possible to verify 
20 that: 

r( r 1) l (60 

This suggests an effective method for the updating of the Cholesky factor 
through the use of unitary transformations, so that, starting from the factor calculated at 
instant n-1 and generator at instant n, it is possible to determine the new factor at 
25 instant n without having to calculate the matrix R{n) . 

In fact, it is possible, arbitrarily assumed L(0) = 0 , to express the updated factor 

as: 

[EW 0] = [E(»-0 g„]Q(") (61) 



L(n)L» = [L(zi-l) g„] 
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where Q{n) is a unitary rotation, in general made as a sequence of elementary 
transformations, Q^Q^A , which, applied to [E(« - l) .g„] produces the nulling of 
the column of the generator g n introducing a zero for each elementary transformation. 

An example of recursive updating, let / 0 (w)and/ 0 (« - l) the first elements in 
position (0,0) of matrices L(n) and L{n - l) , respectively and g n Q the first element of 
the generator vector. It derives that the first line of matrix L(«) can be written as: 

[l 0 (n) 0 0] = [l 0 (n-\) 0 ^ 0 ]Q 0 (62) 
Applying the same transformation to the whole matrix we obtain: 

0 0 



[U«-l) g„]Q 0 = 



'° W L,(»-l) g.. 



(63) 



10 where L,(«-l) and g n[ are the matrices obtained after elimination of the first 
element of the generator column. 

If we recursively proceed in the updating of the whole factor L(n) , we can see it 
factorized in the following triangular form; 

^L(n) 0 



(64) 



_U(«) L»_ 

15 where L and U are the matrices defined in the (53), and L, is the Cholesky factor of 
R^, such for which: L,L* = R^. As it can be noticed from the (64), matrices L and 
L, are triangular while the U is necessarily rectangular, or possibly squared. 

The factorization of the (64) results being essential since it enables the 
calculation of the vector of weights w opt (n) without making the inversion of the matrix 

r (see equation 50) and, as we shall see in short, without inverting the matrix R^ , 
being this operation requested by the algorithm of inverse powers for the calculation of 
the eigenvector ( h opt ) associated to the minimum eigenvalue. 

The recursive factorization for the updating of the whole Cholesky factor requires 
a number of complex multiplication equal to 2(M + L){M + L - 1 ) and a number of sums 
25 equal to (M + L){M + L - 1 ). 
Step 2 

The method of inverse powers for the calculation of the eigenvector associated 
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to the minimum eigenvalue requires the iterative solution of the equation: 

h„ = (R7(n))'\ opl (n-l) % (65) 

h.» - ^ (66) 

n ft 

where the normalization of the (66) is introduced as the result of the constraint on 
power of the impulse response of the estimated channel: |h| 2 = 1. Due to the above 
mentioned constraint, it is convenient to select an initial vector h opt (6) with at least one 
non null element. 

The Cholesky factorization enables the calculation of the above mentioned 
eigenvector without having recourse to the inversion of the matrix, since the (65) can 
be re- written as: 

L>)L>)h B = h op >-l) (67) 
The equation (67) can be factorized in two elementary phases, introducing an 
intermediate vector u„ defined as: 

L»h„ = u„ (68) 
in order to be able to express u rt as solution of: 

L»u„ = h 0 >-l) (69) 
We can solve the equation (69) regarding u„ and then, having found u„ , we can 
solve the equation (68) regarding h n . Both the equations are simple to solve since the 
concerned matrices are triangular matrices. Should we have to solve equations with 
generic square matrices L x L, we should expect a number of arithmetical operation in 
the range of L 3 . On the contrary, being the matrices triangular it is possible to reach the 
solution in only L 2 operations. To illustrate this approach, known in literature and called 
back-substitution, we develop equations (69) and (68): 

l 2x u x +l 22 u 2 = h opt2 (n-\) 
A = A 

l LX u x +l L2 u 2 +h+l Ll u L = h optL {n-l) 
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1 'llK.L = «L 

l LL-A.L +l l-\.L-\ h n.L-\ = "i-l 
A = A 

JlA.L + ClA.L-l +A+/ l"A.I = «l 

Starting from the first term of the (70) we can calculate «, 



(71) 



2/, = 



and rectace it in the other L - 1 equations. This demonstrates that we can reduce the 
problem to a triangular system of dimensions (L - 1) x (L - 1) at the cost of L - 1 
multiplication-and-sum operations and one division. With this method we can solve the 
system at tne cost of L divisions and (L - 1) + (L - 2) + ... + 1 multiplication-and-sum 
operation, for a total of ViL(L - 1) multipiication-and-sum operations and L divisions. 
After having solved the equation (69) regarding u rt we can solve the equation (68). 
The number of operations required is fully identical to that required for the solution of 
the first triangular system of equations. 
Step 3 

Also the calculation of the optimum coefficients of the beamformer can be made 
for a back-substitution in (56), directly in the space of the signals (without explicit 
calculation of the covariance matrix), calculating: 



L »w Q An) = U>)h opf (*) 



(56) 



In this case the previous step 2 must b extended introducing also the calculation 
of the second term of the (56), for a total number of ML + VzU2 multiplication-and-sum 
operations and M divisions. 

The following table summarizes data relevant to the calculation complexity in the 
three steps described, giving also the calculation power of an hypothetical processor 
(expressed in MIPS) requested for the real time execution. It is considered an updating 
of coefficients, extended to the sole midamble, for an array of antennas of M = 8 
receivers and for a matched filter having length L = 5. 



Step no. 


No. of complex multiplications 


No. of complex sums 


MIPS 


1 


2(M+ L)(M + L- 1) 


(M+L)(M + L- 1) 


8,97 


2 


L(L- 1) + 2 L 


UL-1) 


0,86 


3 


ML + YtM 2 


ML + 'AM 2 


0,098 
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The imporovement introduced by a joint estimate can be demonstrated 
considering the previous scenario proposed in fig. 4, where the time dispersion 
introduced by the channel is represented by an echo of the signal delayed by 5 /js. 

Fig. 7 shows the radiation diagram on the azimuth plane determined by the 
beamformer BMF of fig. 6, where it is clear the different behaviour versus the situation 
of fig. 4 concerning the delayed echo. The radiation diagram of fig. 7 shows in fact a 
power level in the echo direction, of the same extent of the level in the useful direction. 
This enables, to the price of a controlled intersymbols interference, to maintain high the 
power of the useful signal present in the estimate process of the transmitted sequence 



The robustness of the new functional proposed can be demonstrated in fig. 8, 
just considering the behaviour of the beamformer BMF in the same critical scenario of 
fig. 5, where a replica of the co-channel interferer was present in the same direction of 
the useful signal. 

As it can be noticed from the radiation diagram of fig. 8, the power is now mainly 
transferred in the useful echo direction, so that the diagram shows a zero also in the 
useful direction, in addition to that of the interferer. This behaviour can be explained 
since the introduction of the time dimension in the functional enables the combination 
of echoes, even when these should be time delayed. We can notice the capability of 
the mixed functional to select in space, among the different delayed echoes that are 
coherent with the training sequence, the one characterized by the greater ratio 
between signal and interferer. 

This behaviour enables to manage also mobile users with superimposed nominal 
directions, whenever they have an opportune spatial dispersion of echoes. It can be 
experimentally demonstrated how, in RCS environment, this enables a considerable 
reduction of the intercell handover probability. 

The joint estimate of the optimum equalization parameters, however requires a 
higher number of parameters than a time only estimate. In particular, in the case of 
FDMA-TDMA systems with comparatively short training sequences Gust as in the 
GSM), two potential problems arise: 

• a growth in calculation complexity in the joint equalization method, compared to a 
disjoint equalization of the prior art; 

• a higher variance in the estimate of the optimum equalization parameters, and 
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therefore a slower convergence speed of the same. 

The length of the training sequence represents a critical parameter in the 
equalizers used, or that can be used, , in the above mentioned systems, since the 
algorithms for the joint estimate of the optimum equalization parameters must identify 
5 the channel impulse response, filtered by the beamformer BMF, availing of a limited 
number of samples, that is the 26 of the midamble of the GSM burst. 

Fig. 9 shows, from top to bottom, a sequence of radiation diagrams on the 
azimuth plane and of representations of the channel impulse response, "photographed" 
on completion of n iterations of their joint estimate, made through the algorithms 
10 adopted at steps 1 , 2 and 3. The scenario shown relates to a useful signal and to a co- 
channei interferer (denoted with "o M and V respectively on the outer arc), both 
accompanied by their own echo delayed by 10 jus (denoted with the same symbols 
placed just under the outer arc). 

Making reference to fig. 9 it is possible to notice how the total convergence is 
15 dominated by the convergence speed of the impulse response h opf («) . At instant n = 

6, when the statistical basis is not sufficiently extended to make a correct estimate of 
the impulse response, the beamformer BMF introduces some zeroes also coinciding 
with the echoes of the useful signal. Also for n = 16, that is when the estimate of 
h {n) leads to recognize the sole main impulse, the beamformer BMF does not 

20 behave yet in an optimum way, since it introduces a zero in a position coinciding wjth 
the delayed echo of the signal. However, we can already notice that for n = 26, as well 
as for n = 30, the channel identification enables the identification of the delayed echo 
as useful signal, and therefore the radiation diagram shows the zeroes only coinciding 
with the interferer and the relevant echo. 

25 A more detailed analysis of the convergence speed can be made evaluating the 

BER at the output of the receiver of fig. 6 as the midamble length varies, when it 
operates in critical conditions, for instance with C/l - - 10 dB. 

Figures 10 and 11 show the results of this analysis referred to a scenario that 
includes a co-channel interferer, in addition to the useful, for three different values of 

30 the signal-noise ratio SNR at the input of the receiver: 10, 20, and 30 dB. More in 
particular, the curves have the length of the training sequence on the x-axis, expressed 
as number of samples, and the BER on the y-axis. Fig. 10 is referred to a situation with 
C/l = - 1 0 dB, while fig. 1 1 to a C/l = 0 dB. 

Making reference to figures 10 and 11, the receiver performance can be 
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considered stabilized in both the cases for values of midamble length between 20 and 
26 samples. In fact, the almost total slope variation of the curves is in this area, and 
this indicates that the convergence of the optimum equalization parameters inside the 
training sequence is reached. 

Making reference to fig. 12, we notice that the general block diagram of a 
receiver whose equalizer is the one shown in fig. 6 with the addition of a portion that 
considers the adaptive operation out of midamble. 

Comparing figures 6 and 12, it can be noticed that the two representations of the 
equalizer differ due to the fact that in fig. 12: 

- and additional RIT block is highlighted, introducing a delay D on the path of the bus 
of the input vector y(n) towards the FUNZ block; 

- the FUNZ block is reached by a vector x!(n) instead of the training sequence {x n } ; 

- there is a new DECTEMP block that stored the vector 3(n), coming out from the 
STIM block. 

15 As we shall see describing the relevant operation, the schematization of fig. 12 is 

suitable to the operation in presence of channels with time varying characteristics 
and/or in presence of non-stationary interferers, that is enabling out of the training 
sequence, which therefore cannot be identified and cancelled considering the sole 
midamble in the estimate of the optimum parameters of the equalizer. Just for this 

20 reason the two schematizations of figures 6 and 12 are complementary between them, 
meaning that during the processing of the midamble the first one is active, and at the 
end, and for the whole remaining duration of the burst the second one is active. In 
practice both the schematizations can be referred to two different operation modes of a 
unique structure of equalizer. 

25 As it is known the Doppler effect is one of the causes that in a mobile radio 

context more contributes to render the characteristics of channels time varying. It 
consists in a small frequency variation of the signal following the movement of the 
mobile, which involves a slow time variation of the amplitude of echoes received. The 
channel thus results time-varying, and even if the variations take place very slowly, can 

30 be useful to try to follow them updating time by time the optimum values of parameters 
h and w estimated during the midamble. In this way we can improve the 

**opt opt w 

demodulation quality, since the estimator STIM is made adaptive, since it recalculates 
at each step the branch metrics on the basis of the estimate of h opt , which in its turn is 
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affected by the w^, one. 

From the operational point of view, it may be convenierrt to update Recording to 
the method described later on) the optimum couple of parameters h opt and vt opt 
estimated during the midamble, and simultaneously equalize the signal received, 

5 starting from the two ends of the midamble and proceeding in parallel in the two 
directions towards the burst ends, previously stored and synchronized. The deriving 
higher circuit parallelism enables to avail of a double time for the execution of the 
calculations. This is a significant advantage, considering the high number of single 
calculation elements to perform in a symbol time of 3,69 us only, which would be 

10 doubled using two semibursts in parallel. 

The presence of non-stationary interferers inside the transmission bus can be 
explained by the fact that the BS relevant to different cells are not necessarily 
synchronized among them. For instance, in the case of an RCS approach to the 
problem of the increase of the spectral efficiency, we can expect to have co-channel 

15 interferers in cells very close to that where is the mobile transmitting the useful signal. 
Since it is not possible to synchronise all the BS, it is possible that an interferer disturbs 
the transmission only during a part of the burst and therefore it cannot be "identified" 
and cancelled during the transmission of the training sequence. 

The equalizer of fig. 12, thanks to the possibility to update during the whole burst 

20 the couple of optimum parameters h opt and w opt estimated during the midamble, is 
thus enabled to face also the case in which one or more interferers disturb the useful 
signal starting from an instant outside the time slot where the training sequence is 
transmitted. 

In the operation, we can say just from now that, with the details we shall give 
25 below, the algorithm developed by the FUNZ block remains almost unchanged. This 
can also be understood from the great similarity existing between the pictures in figures 
6 and 12, made more evident due to the fact that the RIT block, delaying the vector 
y{n) of a discrete time D can be included in the FUNZ block. 

The adaptive operation requires the knowledge of the sequence of data 
30 transmitted {*„} . When these are not known, that is when the training sequence is no 
more available, the same must be replaced by data {B n } decided in the demodulation 
process. It is thus introduced an inevitable delay between the samples of the signals to 
the antennas and the corresponding data of the MLSE estimated sequence, both 
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necessary to determine the error. In fact, as it can be noticed in fig. 12, before the data 
B n relevant to instant n is available, we must wait both for the delay L introduced by 
the matched filter FADAT, and the delay due to the latency time of the Viterbi algorithm „ 
necessary for the emission of sufficiently reliable decisions. More in particular, if D is 
5 the total delay {L + latency of Viterbi), data decided of the sequence {B n } are available 
only at instant n + D , and therefore the estimate of the best equalization parameters 
h (n) and w in) at instant n can avail of the knowledge of inputs on the antenna 

opt v ' opt v ' 

" and of data decided only up to instant n-D, that is the vectors y{n - D) and 
\{n - D) . The RIT block compensates the delay already existing on the other way. 

10 it is worth to say again that vector x(n) contains now the decided data 5„ , while 

during the combined estimate of the optimum equalization parameters, made during 
the midambie, the same vector contained already known data x n and therefore without 
delay. The statement made on vector y{n - D) means that, concerning the calculation 
of h op ,{n) and w op ,(«) , the mean at the n -th instant calculated on the N elements of 

15 the columns of matrices R JB ,R JV ,an</R XK concerns elements included between the 

instant n- D and instant n- {D + N) . 

What in pratice occurs due to the delay D , is that the combined estimate of the 
optimum equalization parameters at instant n results in delay compared to the channel 
variations, which have effect without any delay on symbols y(n) received. Two needs 

20 therefore conflict. The first one is to maintain high the latency of the Viterbi algorithm, to 
assure correct decisions avoiding error propagation in the estimate loop. The second 
need is to furnish the data decided in the shortest time to the FUNZ block, in order not 
to introduce an excessive delay in the updating of the best equalization parameters, 
such to compromise the possibility to follow quick variations in the channel 

25 characteristics. 

The approach followed in the subject method is to use in the estimate of the 
optimum equalization parameters, a. sequence of low delay estimated data, 
represented by the L elements of the vector 3{n). The latter, is obtained forcing the 
estimator STIM to decide just after few steps ( 4 or 5 approx.), which is the more 

30 probable sequence, obtained climbing the trellis the maximum likelihood path. Data 
finally decided by the sequence {^tj)} are on the contrary available at the end of the 
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burst, in order the same are as reliable as possible. 

The above mentioned technique, which solves the problem to avail in advance of 
the data of an almost correctly estimated sequence, .is called "Per-Survivor 
Processing", or PSP technique, and is described in the paper under the title "Per- 
5 survivor processing: a general approach to MLSE in uncertain environments", authors: 
R. Raheli, A. Polydoros and C. K. Tzou, published in IEEE Trans. Commun., vol. 43, p. 
354-364, February/March/April 1995. With this technique the Viterbi algorithm can be 
used in presence of one or more unknown, or not well determined parameters, in this 
case the samples of the optimized impulse response, which if on the contrary were 

10 known could be profitably used in the updating of branch metrics. 

Since low delay decisions can result wrong, an error propagation phenomenon 
can be originated due to a vicious circle establishing in the loop formed by the 
feedback of 3(n) toward the block FUNZ. In fact, considering all the elements of 
vector S(n) as the result of single low delay decisions at symbol time, as belonging to 

15 a sequence running in a shift register long L (this is what actually occurs in the 
representation of the states according to Viterbi), .a wrong decision would finally 
produce its negative effects for a number of symbol times equal to that defined for low 
delay decisions. 

The subject equalization method eliminates also this second drawback by 
20 replacing, at any symbol time «, all the L elements of vector j!(n) by a same number 
of temporary data . directly copied from the maximum likelihood path in the trellis valid 
at the time of n-th symbol. 

Fig. 13 shows in fact an example where at instant N a wrong temporary decision 
is taken in the trellis, which however does not propagate because at the subsequent 
25 instant N + 1 the decided sequence copied is the correct one. Actually, not always 
things can proceed like this, but very likely this is verified in presence of interferes not 
particularly intense. 

The operation method just described makes the updating of the optimum 
parameters out of the midamble more accurate, remaining the approximation 
30 introduced by low delay decisions. Furthermore, to improve the adaptive behaviour of 
the equalizer of fig. 12, it is necessary to introduce a X coefficient so-called "forgetting 
factor", in order to less weigh the/ess recent data. By this, recalling the (58), the matrix 
R at instant n can be written as: 
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R(n) = m(n-\)~g n g n (72) 

"yOO 



15 



o = 



=w ^ (73) 

where A is an arbitrary real constant for which: 0 < A < 1 . and g„ is an updating vector 
defined as in the (59), with the difference that its elements are here the low delay 
decisions. 

Likewise, recalling the difference matrix (57) and the (58), we can write: 

R(n)= F(«)R(n - l)F*(n) + g„g'„ (74) 

being: 

F(«) = VTl (75) 
With this position the updating of the Cholesky factor can be directly obtained 
from equation (74). In fact, let L the Cholesky factor of R , we can update the same, 
likewise we have done in the (61), according to the relation: 

[L(«) 0] = [VAL(zi-l) g„]Q(«) (76) 
With all the above clarifications, the estimate of parameters h opt {n) and w opf (n), 
made adaptive, is obtained applying the recursive algorithm already described for the 
non adaptive estimate made by the FUNZ block of fig. 6. 

Making reference to figures 14, 15, 16, and 17, the actual adaptive equalization 
of the subject method are shown. 

Fig. 14 shows two BER trends along the burst of a same number of demodulated 
20 signals, for an interferer disturbing up to the 50th symbol. A first trend (dashed line) is 
relevant to the subject equalisation method, assuming to 

freeze the estimate out of the midamble, while the second trend (solid line) relates to 
the same method enabling the performance of the adaptive operation. The input signal 

E C 

snows ±±_ _ iQ & anc j — = 0 dB ; the trends are averaged on 500 performances. The 

Wo I 

25 BER indicates the probability that the generic symbol is demodulated in a wrong way. 
for instance for bad equalization. As it can be noticed from the first curve, if the burst is 
demodulated using the equalization parameters estimated with the training sequence, 
the performance greatly deteriorate for the time the interferer is present. 

• The second curve shows on the contrary, that the adaptive operation can limit 

30 the degradation of the performances. In fact, proceeding from centre to the ends, it can 
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be noticed how the performances deteriorate coinciding with the arrival (in the 
calculation) of the interferer, to better later on, after equalization. This means that the 
equalizer could adapt to the new transmission situation and vary the radiation diagram 
of the antenna array in order to cancel the interferes The next figures 15, 16, and 17 
5 show radiation diagrams on the azimuth plane relevant to the curve with solid line of 
fig. 14. 

Fig. 15 shows the radiation diagram estimated at the end of the midamble, when 
the interferer is not still present. The arrival directions of the echoes of the useful signal 
are represented with solid line with line lengths having length proportional to the power 
10 of the relevant echo, expressed in dB. 

Fig. 16 shows the evolution of the previous radiation diagram immediately after 
the identification of the interferer, whose echoes are represented by dashed lines 
having length proportional to the power of the relevant echo, expressed in dB. 

Fig. 17 shows the radiation diagram when the joint estimate of the optimum 
15 equalization parameters has reached the convergence. As it can be noticed, the 
beamformer can cancel the signal coinciding with the arrival directions of the 
interferers. 

The evaluation of figures 14 and 17 shows that applying the equalization method 
scope of the present invention, the joint estimate of the optimum equalization 
20 parameters can converge towards a couple of vectors h opl {n) and w op/ (w) that can 

assure a good quality in demodulation, even in presence of non-stationary interferers. 
What just stated is true until the non-stationary interferer is in the same size range of 
the useful signal. In the contrary instance, a sudden peak of the demodulated signal 
can occur in the BER, that can hardly be taken to acceptable values in a reasonably 
25 short time compared to the residual duration of the burst. 

Fig. 18 shows three trends of BER along the burst of a same number of 
demodulated signals, for an interferer disturbing up to the 50th symbol. The three 

E b 

curves are referred to an input signals with T}~= 10 dB and three different values of 

C 

— , of 0, - 10, -20 dB, respectively; the trends are averaged over 500 performances. 

30 The worsening of the convergence as the ratio C/I decreases is evident, in particular, 
we see that the BER does not drop at acceptable values until the ratio C/I remains 
lower than -10 dB. 
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The cause for this behaviour can be due to the delay, even small, affecting the 
low delay decisions in the MLSE estimate of data transmitted. * 

Just due to this delay the equalizer does not immediately detect the arrival of a 
strong interferer, in such a way that, when the first temporary decisions that consider 
5 the interferer are available, therefore suitable to take measure to cancel the same, it 
will have in the meantime caused a great number of wrong decisions, having an extent 
such to compromise the convergence of the estimate made by the FUNZ block. 

The inconvenient is highly attenuated by a method alternative embodiment, 
which tries to cancel the delayed vision of the interferer in the adaptive equalization, in 
10 order to quickly cancel the non-stationary strong interferes. 

The embodiment includes a phase for the continuous monitoring of power 
received on the antennas, and for the detection of possible peaks higher than a set 
threshold. The detection of a peak at an instant n is used to automatically select the 
operation of the equalizer according to the methods of the embodiment. Since noises 
15 to be detected are sufficiently strong ones only, a particularly sophisticated power peak 
detection method is not necessary. For instance, it is possible to use a non-coherent 
detection method, known to those skilled in the art. In the receiver of fig. 12, for ease of 
representation, we have not included an appropriate block performing the above 
mentioned functions. 

20 The detection phase is profitably made in parallel to normal operation out of the 

midamble. The embodiment quickly exhausts its effects in a transient time having AT 
duration (of about twenty samples) necessary, as we shall see later on, to the setting in 
the computation of some matrices used by the FUNZ block. When AT elapses, it is 
convenient to automatically reselect the original operation that results more effective at 

25 transient exhausted. 

The strategy implemented by the embodiment to cope with strong intensity 
stationary interferers, is to advance at instant n the future knowledge about the 
interferer to cancel. In this way the interferer is "seen" in advance in the calculations, 
that is before the first relevant decisions are available and the adaptive algorithm of 

30 error minimization can start updating the beamformer coefficients to try to eliminate the 
noise. The above mentioned strategy is implemented in transient AT , dividing in two 
elementary phases (1. and 2.) the updating of matrix R(n) . In practice, we proceed as 
follows: 

1. during this first phase, which elapses within the symbol time n , it is anticipated the 
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updating of the sole spatial covariance matrix R v> , with vectors y{n) obtained from 
the samples of signals on the antennas, present starting.from instant n up to instant 
n + AT in the stored burst. The statistical basis of the above mentioned matrix is 
thus extended, that now contains sufficient information of the interferer to cancel. 
5 The matrix R is no more modified up to the end of the transient. Then we proceed 

to phase two. 

2. During this second phase, for n < n < n + AT , the covariance matrices on R ^ 
data with vector xP(«), and the cross-covariance matrices R ^, with vectors £'(«) 

and y(«) downstream the RIT block are updated. Concerning vector H(n) it is 
10 worthwhile pointing out that starting from « , once the low delay decisions term is 
elapsed and in lack of appropriate measures which shall be shortly specified, this 
vector could express completely wrong decisions due to the performance of phase 
1. 

When the n + AT instant is reached the estimate of the mentioned matrices 
15 proceeds in the conventional manner. The updating process of matrix R(w) , in the time 

interval [« , n + AT\ , can therefore be expressed as follows: 

expressing the different time interval used for R^ calculation with respect to the other 
terms. 

20 The "anomalous" updating of matrix R(n) generates a setting transient having 

duration AT in the computation of matrices R^ and R^ . Said interval just coincides 
with the amplitude of the time shifting in the updating of the spatial covariance matrix 

The joint estimate of the optimum equalization parameters will therefore result 
25 "vitiated" by this type of approach, which however proves to be necessary to cope with 
quick variations of the channel characteristics, enabling the beamformer BMF to 
introduce some zeroes coinciding with the echoes of the new interferer, which 
otherwise it will not be able to introduce, unless with an excessive delay . 

As a consequence of the matrix setting time R(n) too many errors could 



R(«) = 



(77) 
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propagate in the estimate loop during the transient [«, n + AT] , therefore a further 
feature of the method embodiment, is to update the vector w opt {n) of the weights of 
the beamformer BMF during the transient only, maintaining the vector of the samples of 
the impulse response constant for the whole transient time at a value h opt it had at 
instant n . The choice rewards the coefficients of the beamformer BMF since they 
result less sensitive to decision errors compared to the samples of the impulse 
response. 

For what said above, we can now estimate even during the transient 
[«, n + AT] the optimum equalization parameters. In particular, h opt is constant and 
the (51), recalling the (77), enables to determine the vector w^, , solution of the 
equation: 

R w {n + AT)* op , («) = R * 0 , (n)h opl . (78) 

in which R (n + AT) is calculated at instant n and then kept constant for the whole 

AT, while R'Jn) is calculated at each instant «-th" of the transient. As it can be 

noticed, due to h op , constant, the (67) and (68) need no more to be solved. 

The previous considerations on the updating methods of matrix R(n) of the (77), 

in the time interval [n, ii+AT], are verbatim transferred to the updating of the 

Cholesky factor of the same. 

Concerning this last point, at an instant n of the detection of the power peak 

received on the antennas, the following actions are made: 

- within the symbol time n , it is updated the sole triangular matrix L with vectors 
y(n) obtained from the samples of signals on the antennas, present in the stored 
burst, starting from instant n up to instant n + AT ; 

- for n < n < n + AT , matrices U(«) and L» rectangular and triangular, 
respectively, are updated with said updating vector g„ and equation (56) is solved, 
keeping the vector of samples of the impulse response h opl constant at the value 

assumed at said instant n ; 

- the calculation methods before said instant n are restored when time n +AT, 
elapses. 
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L(«) = 



(79) 



For what said above, the Cholesky factor can be submitted to triangular 
factorization in this time interval, as follows: 

'L(n-hAT) 0 
UU) L»_ 

With all the above considerations relevant to the possibility to use in the 

5 calculations either directly the matrix R(«) f or preferably the Cholesky factor L(n) of 

the same, the different steps of the adaptive estimate of the sole vector w opt {n) of 

weights characterizing the method embodiment, are connected to the corresponding 
ones of the adaptive algorithm as for the utilization method in the computation of low 

delay decided data {j3 n } . Said adaptive algorithm in its-turn, referred to the detailed 
10 description of the algorithm developed by the FUNZ block included in the receiver 
sketched in fig. 6, this description continuing to be valid making the modifications 
relevant to the time progress AT of matrix R xy (« + AT) or of the relevant Cholesky 

factor L(>7 + AT) . 

This assumed, the weight vector w op/ («) can be calculated either solving the 

15 (78) or, preferably, through Cholesky factorization and subsequent solution of the (56) 
which in this new context becomes: 

V {n + AT)v, opt {n) = U»h opf (56") 
Concerning this second opportunity however, it is necessary to analyze how to 
update the Cholesky factor, when it assumes the expression of the (79), applying the 
20 unitary transformation Q(w) defined in the (76). In the case of use of the embodiment, 
contrarily to the conditions that applied to the (61) and (76), it is necessary to update, 
at instant n , only the L factor of the spatial covariance matrix R^. Analyzing the 

mechanism according to which the unitary transformations are applied, we can realize 
how it is possible to update matrix L only "delaying" the updating of the other factors. It 

25 was previously said, making reference to (61) that a unitary rotation is generally 
implemented as a sequence of elementary transformations, Q 0 ,Q, A , usually Givens 
rotations. The updating delay can be introduced remembering what expressed in the 
(62) and (63), where the single elementary transformation Q, acts on the single 
column I,(w- l) of factor L and on the column of generator g n . If we apply the. 

30 transformation Q, to the first M elements of the Mh column, where M is the order of 
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matrix L , it is possible to delay the updating of factors U and L J( storing the 
transformations occurred in L and restarting the updating of U and once low 
delay decisions are newly reliable, that is at the end of the setting period AT of these 
matrices. 

5 Fig. 19 shows two trends of BER along the burst of a same number of 

demodulated signals, for a strong interferer disturbing up to the 50th symbol. A first 
trend (dashed line) is relevant to the out-of-midamble adaptive method of the subject 
equalization method, while the second trend (solid line) relates to its alternative 

E 

embodiment for strong interferes. The input signal shows —f-=\QdB and 

10 y = -20 dB ; the trends are averaged over 500 bursts. It can be clearly noticed that the 

new method starts the adaptation first, thus succeeding in containing the peak in the 
trend of the error probability, and assuring a significant improvement in the quality of 
the demodulated signal. 

Further bibliographic references useful to the comprehension of the 
15 "Displacement Theory" applied to the time-variant channel are reported in the following 
papers: 

« "Time-Variant Displacement Structure and Triangular Arrays", by: A. H. Sayed, H. 
Lev-Ari and T. Kailath, published on IEEE Trans, on Signal Processing, vol. 42, No. 
5, May 1994. 

20 • "Fast Computation of Channel-Estimate Based Equalizers in Packet Data 
Transmission", authors: N. AI-Dhahir, J. M. Cioffi, published on IEEE Trans, on 
Signal Processing, vol. 43, No. 11, November 1995. 
• "A Computationally Efficient FIR MMSE-DFE for CCI-lmpaired Dispersive 
Channels", by N. AI-Dhahir, published in IEEE Trans, on Signal Processing, vol. 45, 
25 No. 1, January 1997. 

The teachings of the equalization method scope of the present invention, and its 
alternative embodiment, find the relevant application in the equalizers included in the 
functional blocks of figures 6 and 12. The implementation of the above mentioned 
blocks is known to those skilled in the art, as well as the implementation of the 
30 functions indicated in the text which do not find corresponding blocks in the figures, 
except for the FUNZ block and the relevant internal blocks PESI and CAMP. Moreover, 
as a consequence of the different interrelation among the blocks that a joint estimate of 
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parameters w ,h imposes versus a separate estimate, also the equalizer 

architecture does not result known. K 

In the detailed description of the FUNZ block below, reference shall be made to 
the sole representation of fig. 12, which enables to implement a sole equalizer capable 
5 to operate both during the midambie, and out of it, in this last case according to the 
adaptive method, or according to that of the alternative embodiment, after appropriate 
hardware/software selections controlled by events such as the end of the midambie, or 
the sudden arrival of a strong interferes 

As far as the circuit realization of the FUNZ block is concerned, it is preferred to 

10 have recourse to circuit structures known as "systolic arrays", represented in figures 
21, and 22. These structures have a high parallelism degree, that makes them suitable 
in themselves to operate on the matrices, but mainly to make almost simultaneously a 
number of calculations which should otherwise be made in sequence. A 
microprocessor based implementation would be well less effective, just because a 

15 common DSP is governed by a time sequence of elementary instructions making it 
more suitable to solve serial problems, contrarily to the systolic array that performs all 
the calculations in parallel in one symbol time of 3,69 (is (or in a double time in case of 
joint operation on the semibursts). Even if one wants to have recourse to a DSP of the 
RISC type, the latter would not have a parallelism comparable with that of the systolic 

20 array, which seems to be therefore "problem oriented". In the practice, the FUNZ block, 
as well as the remaining blocks of the equalizer of fig. 12, are profitably implemented 
through digital integrated circuits of the ASIC type (Application Specific Integrated 
Circuit). 

Making reference to fig. 20, it can be noticed that the FUNZ estimate block 
25 includes five differently configured systolic arrays, indicated with SYSTOL1, SYSTOL2, 
SYSTOL3, SYSTOL4, and SYSTOL5, respectively. The vector x(n) referred to the 

samples of the training sequence {x n } stored in the receiver reaches a first input of the 
SYSTOL1 block, a second input is reached by vector y(n) , and the third input by a 
MID control signal resulting true only coinciding with the midambie. On a first output 
30 bus the Cholesky factor L(w) indicated in (61) is present, and reaches a first input of a 
two-input selector SEL, while a second output bus CONT connects SYSTOL1 to a first 
input of SYSTOL2. 

The second input of the SYSTOL2 block is reached by vector 3{n) referred to 
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the data of the low delay decided sequence {B n }\ a third input is reached by vector 
y(n) coming from the RIT delay block; a fourth input BURST is a bus having reading 
access to a memory in which the receives writes the samples y{n) of the reception 
burst; a fifth input is reached by the forgetting factor A , a sixth input is reached by a 
PIC signal, we shall describe below, and finally a seventh input is reached by the MID 
control signal. The Cholesky factor L'{n) is present on the output bus as indicated in 
the (76). that reaches the second input of selector SEL, whose selection input is 
reached by the MID control signal. The bus coming out from selector SEL is connected 
to one end of the switch block SWT controlled by a PIC signal having duration AT 
coming from a PICSEL block. This last block receives the MID signal and two other 
signals indicated with PIC and TRANS on three relevant inputs, and on this basis 
generates the PIC signal. The TRANS signal is active for a time AT equal to the 
duration of the transient operation foreseen for the method embodiment; it is a 
configurable parameter of the equalizer, like the forgetting factor X . The presence of a 
logic, taking care of the generation and control of all the indicated signals, is obvious. 

The PIC signal comes from a PICRIV block receiving the high frequency receipt 
signal directly from the antenna array ARY. The bus coming out from the other end of 
the SWT block is connected to an input of SYSTOL3, from which the intermediate 
vector u„ comes out, directed to a first input of SYSTOL4, whose other input is 
connected to the input of SYSTOL3. A vector h op ,(«- 1) is present at a first output of 
SYSTOL4, directed towards an additional input of SYSTOL3. The vector h^n) is 
present at a second output of SYSTOL4, directed towards a first input of SYSTOL5, 
whose second input is connected to the output of the SEL block. At the output of 
SYSTOL5 the vector w opl {n) is present. The two blocks SYSTOL3 and SYSTOL4 are 
included in a functional block AUTOVET, whose input coincides with that of SYSTOL3 
and whose output coincides with that of SYSTOL4. 

The architecture of the FUNZ block of fig. 20 enables a full operation of the 
equalizer of fig. 12 both during the midamble and out of the same, according to the 
method scope of the invention and its alternative embodiment. 

During the midamble, the MID signal results true, this condition is used to enable 
the array SYSTOL1 and disable SYSTOL2, to control the selector SEL in the selection, 
of L(«) , and determine the closing of the switch SWT. Consequently, the bus of 
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SYSTOL1 is extended towards the AUTOVET block extracting the portion relevant to 
matrix L , and uses it to solve the (68) and (67), determining the vector h opt [n) . More in 
particular, the array SYSTOL3 calculates the intermediate vector u„ according to the 
(68) applying the "back-substitution" method, while the SYSTOL4 array calculates 

5 h (w) according to the (67), again through "back-substitution". The array SYSTOL5 
extracts from the bus coming from SEL the portion relevant to matrices handV and 
uses it, together with vector h 0/w (w), to solve the (56) and determine the 
vector w uw (/7) , always via "back-substitution". 

During the adaptive operation, the negate MID signal disables the array 

10 SYSTOL1 and enables SYSTOL2, and controls the selector SEL to the choice of 
L'(rt). In lack of a change in the logic condition in the PICsignal, the switch SWT 
remains closed and the bus of SYSTOL2 is extended towards the AUTOVET and 
SYSTOL5 blocks; these last treat the relevant matrix components of the Cholesky 
factor L'(«) like what said for L(«). 

15 Before the SYSTOL2 array starts operating, the transition from true to negate of 

the MID signal controls the transfer of the content of array SYSTOL1 on the CONT bus 
to be copied in the cells of SYSTOL2. The above mentioned operation enables the 
updating of the already existing Cholesky factor L(w) according to the new calculation 
method foreseeing the use of the forgetting factor X. If the above transfer were not 

20 foreseen, the array SYSTOL2 would newly estimate the factor L'(w) but with lower 
possibilities versus the previous situation. 

When the PICRIV block detects, at an instant n , the presence of a power peak 
on the antennas, higher than the set threshold, the PICSEL block discriminates 
between two possible operation conditions. A first condition is that of strong interferer 

25 in presence of midamble, in this case the PIC signal results false and the estimate 
during the midamble is not disturbed. The second condition is the strong interferer one 
out of the midamble, in this case the true signal PIC 1 controls the opening of the switch 
SWT, and controls appropriate switches included in the array SYSTOL2 that make a 
temporary modification in the structure of the array such to render it suitable to operate 

30 as foreseen in the method alternative. With true PIC the AUTOVET block does not 
receive anything at input and freezes the value of vector h op , at instant n up to 
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elapsing of A7\ The SYSTOL5 block solves the (56), always through 
"back-substitution", using the vector h opl and the matrices LandU of the Cholesky 

factor L'(n) derived from the above mentioned structural modifications made to the 
array SYSTOL2, determining in this way the vector w op /(w) . 

5 When time n -f elapses, the new condition of the PIC signal controls the 

closing of the switch SWT, but the configuration of the array SYSTOL2 remains as 
modified up to the end of the burst, and the equalizer operation returns to the adaptive 
method. 

Concerning the circuitry realization of the different blocks of fig. 20, the following 
10 must be pointed out: 

- blocks PICRIV, PICSEL, RIT, SEL, and SWT, on the basis of the functional 
description supplied, can be implemented applying the common technical 
knowledge; 

- the systolic arrays SYSTOL1, SYSTOL3, SYSTOL4, and SYSTOL5 can be 
15 represented by an architecture of cells common to that of all known systolic arrays 

of triangular form. Their description therefore requires to specify the nature of the 
cells, to highlight where the array diversify among them and from those of the prior 
art. A functional description of the cells indicating the type of arithmetical operations 
made on the entering signal and on the pre-existing content to obtain the new 
20 content and the output signal, must be considered sufficient to their realization 
availing of the common technical knowledge in the field of systolic arrays and gate 
arrays. 

- the systolic array SYSTOL2 avail of the same structure of SYSTOL1, but includes 
FIFO memories and electronic switches in the descent chain, as we shall see in fig. 

25 22. The ways to obtain the above mentioned FIFOs and relevant switches in chips, 
are still those that enable to obtain the different arrays, that is the gate array 
technology of the ASIC type. 

Of course, the FUNZ block can be implemented through architectures differing 
from the indicated one, it is possible for instance to use a unique systolic array 
30 grouping the SYSTOLI and SYSTOL2 functions, in this case it should be useless the 
transfer of the content to the end of the midambie. This new array must however be 
able to modify at the end of the midambie the arithmetic operations defining the content 
of the cells, since now it must be considered the forgetting factor X . 

With reference to fig. 21 , in which the systolic array SYSTOL1 is represented in 
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functional blocks, we can notice that there are single blocks having square shape, 
interconnected among them, denoted with 1. , ii i e I , and peripheral bjocks having 
circular shape c } connected to the first ones. The structure of the array has triangular 
shape, as that of the matrix L(w) it reproduces. Actually, each block J i , u i , I s { , and c i 

5 represents an element of this matrix. As it can be noticed, blocks J iy c i are grouped in 
the area of the upper sharp vertex corresponding to the position of the triangular matrix 
L , blocks l s , c ; are grouped in the area of the lower sharp vertex corresponding to 

the position of the triangular matrix L^., and finally blocks u i are grouped in the 
remaining square area (or, more generally, rectangular) corresponding to the position 

10 of matrix U . Circular blocks c i form the closing diagonal of the array, and the two 
blocks at the ends indicate the starting and the end of the updating. For ease of 
representation 3x3 matrices are shown. 

Fig. 22 indicates the architecture of array SYSTOL2, obtained from that of fig. 21 
interrupting the output line deriving from the blocks corresponding to the elements of 

15 the last line of the triangular matrix L , including identical FIFO memories having length 
such to introduce a delay equal to the transient AT , and connecting the outputs of the 
FIFOs to relevant blocks corresponding to the elements of the first line of matrix U . 
The array SYSTOL2 can operate according to the adaptive method and to the method 
embodiment, the means enabling this possibility shall be described. 

20 Fig, 23 shows an electronic switch COM, identical to the other ones included in 

the array SYSTOL2, inserted between two points A and B in fig. 22, corresponding to 
the input and output ports of a generic FIFO. The central point of the switch COM is 
connected to point A of origin of the signal to switch, while the FIFO is connected 
between the two output terminals of COM, one of which is connected to point B. The 

25 active PIC signal controls the switching in order to insert the FIFO between A and B as 
shown by the position of fig. 23; the signal propagating on the corresponding column is 
thus delayed by AT . However, at the end of the transient AT , the FIFO remains 
inserted up to the end of the burst, instant at which occurs the opposite switching that 
shortcircuits points A and B, and the delay is excluded returning to the adaptive 

30 operation in presence of a new burst. In the figure it is also highlighted a clock signal 
CKT at symbol frequency that sets the operation time of the FIFO. There will be as 
many COM blocks as are the columns over the FIFOs. 
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Referring to fig. 24, we can notice that circular blocks are affected by an input 
signal in and by an output signal p. Also square blocks are v affected by the same 
signals, but they have also an additional input signal p and an additional output signal 
out. An exception are the square blocks of the last horizontal row and the last circular 
block of the same, where the output signal p is missing. The content of blocks, 
similarly for each type of block, is denoted with a. 

SYSTOL1 and SYSTOL2 arrays are arranged in the figures in such a way that 
input signals, corresponding to the updating vector: 

y(n) j 



<T — 



(n)\ 

io enter from the left side. More in particular, numeric signals y l (n),...,y i (n) ,y,,(«) 

are directed towards the blocks of the first vertical line of the array, corresponding to 
the elements of the first column of the triangular matrix L , while numeric signals 
(«),..., Xf (n),....,x L (n) are directed towards the blocks corresponding to the 
elements of the first column of matrix (rectangular) U . Concerning the SYSTOL2 
15. array, the vector x(«) both during the adaptive operation which during the alternative 
coincides with vector vector y{n) during the adaptive operation shall be that 

coming out from the RIT block, that is y{n - D) ; during the variant on the contrary it 
shall be a vector y{n+AT) consisting of samples 
y,{n + AT),...,y,(n + AT),....,y M {n + AT) of the stored burst, transferred in 
SYSTOL2 and there used in the advanced updating of the Cholesky factor 
L(« + AT) of the spatial covariance matrix R w (b + AF) . 

Even if not expressly indicated, provisions are made for a general clock signal of 
the array of figures 21 and 22, at a frequency much higher than the symbol one, which 
simultaneously reaches all the blocks l t , w, , /,,. , and c, to control the relevant single 

25 operation, the bench of FIFO memories of fig. 22 is updated at symbol time. 

The geometric arrangement of the blocks shown in figures 21 and 22 is observed 
on the chip area, where the blocks form a same number of cells; this enables the signal 
originating in the upper vertex to propagate along the descent chain. It is obvious that 
the correct content of the cells shall be available only after a time sufficient to the 

30 complete propagation of the signal towards the opposite end of the diagonal made of 
the circular blocks, following the longer path between the cells. The single signals 
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entering the left side of the array are associated to a corresponding number of arrows 
that indicates the delay according to which their effect develops. 

In the operation, systolic arrays SYSTOL1 and SYSTOL2 are respectively 
referred to the updating of the Cholesky factor L(w) of matrix R(«) as indicated in the 
(61 ) and (76), we repropose for ease of reading: 

[L(*) 0] = [l(*-1) g„]Q(n) (61) 

[E(/i) 0] = [VAL(«-l) g fl ]Q(«) (76) 
During the strong interferes the updating indicated by the (76) is modified by 
the introduction of the FIFO as seen in figures 22 and 23, but the relations we shall 
now indicate will continue to have value on the elements of matrix L(« + AT) . 

The unitary rotation Q(«) is made associating to each cell a relevant operation 
made on input signals and on the already existing content, to obtain the output signals 
and the new content. In particular, for instance, of the (76): 

- in signals propagate in horizontal direction and correspond to the 
15 jv, ("),..., X-M yM e x,(«),-.., *,(«),..■ signals at the input of the 

cells of the first column of the array; for the remaining cells indicated by square 
blocks the in signals correspond to the out signals coming out from the previous 
cells in the horizontal propagation direction of the signal; 

- signals p propagate in vertical direction and are generated at the output of each 
20 cell, both square and circular, calculating the following expression: 



10 



25 



p = -^= (80) 

where X is the forgetting factor, and a is the new content of the cells calculated as 
follows, concerning the sole circular cells corresponding to the main diagonal of 

matrix L : 



a <_ Jx\ a \> +\in\ 2 (81) 
the new a being pointed by the arrow; 

- for the remaining square cells the new a is determined jointly to the signal out that 
propagates in horizontal direction, in fact, the following matrix relation applies: 



[a out\ <— [aVX j'w] i — 



p -i 



(82) 
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where the vector pointed by the arrow indicates the new values of a e out . 

Relations similar to the (80), (81), and (82) are valid for the updating according to 

the (61), giving A = 1. 

Concerning fig. 22, the introduction of FIFO memories in the descent chain of the 

5 systolic array is the simple circuit translation of what said describing the method 
embodiment concerning the updating of the sole Cholesky factor L(n + AT) . In fact, 
we said that it is possible to delay the updating of U and L s factors storing the results 
of the elementary transformations Q.. The delay introduced at instant n is then 
recovered at instant n + AT , that is when the estimator STIM is able to supply the 

10 reliable symbols decided and can therefore restart a normal updating of the Cholesky 

factor L{n) . 

As for the functional description of the cells belonging to the systolic arrays 
SYSTOL3, SYSTOL4 and SYSTOL5, that is the mathematical relations developed in 
the cells on the basis of the previous content and of new inputs to generate new 

15 outputs, the sector technician is enabled to obtain the above mentioned systolic arrays, 
applying the steps of the back-substitution method described in the text to solve the 
relevant equations indicated, and availing of the common technical knowledge in the 
field of systolic arrays. 

Profitably, the equalizer belonging to the receiver of fig. 12 (in practice the whole 

20 block structure shown in the figure) is duplicated to process in parallel the two 
semibursts at two sides of the midamble. A burst reconstruction circuit placed 
downstream the two equalizers puts in the right time order the data of the two 
estimated subsequences, obtaining a unique sequence {4} of data transmitted 
estimated at maximum likelihood. 
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CLAIMS 

1. Equalization method for the cancellation of isofrequential interferes in 
signals received by array antennas (ARY) of base stations of a cellular mobile system, 
said signals (y{«)) being received under the form of bursts of symbols of a 

5 transmission burst, including a training sequence ({*„}) known to the receiver, said 
bursts received being demodulated, digitized, stored, and synchronized through 
discrimination of the training sequence, characterized in that it includes the following 
steps: 

a) spatial filtering of signals coming from the antennas (y(«)) , using a spatial filter 
10 (BMF) combining among them on the basis of appropriate weights (w opf («)) the 

corresponding symbols in the different bursts stored; 

b) time filtering of the signal {y oul (n)) generated by said spatial filter (BMF) using a 
matched filter (FADAT) by samples (h o;w (w)) of the impulse response of the 
transmission channel; 

15 c) joint optimization of said weights (vf opl {n) y h opt {nj) on the basis of the symbols of 

said training sequence ({*„}) and of said signals coming from the antennas (y(«)) 
present in stored bursts; 
d) estimate of symbols {{B n }) of the transmission burst corresponding to the sequence 
coming out from the time matched filter (FADAT); 
20 e) updating of optimized values of said weights (w op ,(*), h 0/3 ,(«)) . out of the time slot 
of the training sequence, through repetition of the previous steps a)...d), replacing in 
step c) the symbols of the training sequence ({*„}) by symbols estimated ({3 rt }) 

in step d), thus rendering said equalization adaptive. 

2. Equalization method according to claim 1, characterized in that said joint 

25 optimization of weights (w^w), h opt {n)) is made minimizing a functional (FUNZ) 
corresponding to the mean square error originating from the comparison between the 
corresponding samples *-th of two sequences, out of which, a first one is obtained 
submitting the symbols of the training sequence ({*„}) and a time filter (CAMP) that 
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estimates said samples of the impulse response of the transmission channel (h(«)) , 
and a second sequence is obtained submitting the «-th samples of said stored bursts 
(y{n)) to a spatial filter (PESI) estimating said weights {w{n)) . 

3. Equalization method according to claim 1 or 2, characterized in that said 
5 joint optimization of weights converge in a time slot of said training sequence 

having length lower than or equal to 26 symbol intervals. 

4. Equalization method according to any of the previous claims, characterized 

in that when said training sequence ({*„}) is included in a midamble of the 
transmission burst, said equalization starts from the ends of the midamble and 
10 proceeds in parallel in the two directions of the burst, previously stored and 
synchronized. 

5. Equalization method according to any claim 1 to 3, characterized in that 
when said training sequence ({*„}) is included in a midamble of the transmission 
burst, said equalization starts from a first end of the midamble and proceeds towards 

15 the end of the burst, previously stored and synchronized, and afterwards it restarts 
from the second end of the midamble towards the other end of the burst, using 
estimated data of the first semiburst to have a longer group of estimated data for the 
equalization of the second semiburst. 

6. Equalization method according to any claim 2 to 5, characterized in that 
20 said joint optimization of weights and samples made in said step c) includes the 

following steps: 

d1) assignment of arbitrary values, preferably null, to the elements of a structured 
triangular matrix L(«), defines as follows: L(«) = 

L(«) and L,{n) are triangular matrices and U(n) is a square matrix, or more in 

25 general a rectangular one; 

d2) assignment of arbitrary values, one of which non null, to the elements of a vector 
h op , («) storing said optimized samples of the impulse response of the transmission 

channel; 

d3) assignment of arbitrary values, one of which at least non null, to the elements of 
30 a vector w„ ,(«) storing said weights of the spatial filter (BMF); 



L(«) 
U(») 



0 



wherein 
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d4) updating of said structured matrix L(w) applying at instant n -th the following 
iterative relation: [E(w) 0] = - 1) g„]Q(«), where g„ is the following 



updating vector g n = 



y(»)' 
%(«) 



structured in two component vectors, out of which: 



y{n) = [y ] (n\y 2 (n),....,y M (n)Y is a first vector whose M elements are symbols 
coming from M antennas at instant n -th t present in a same number of stored bursts, 
and x{n) = [xS n )> x r,-i( n )>''--> x n+L-\( n )] T is a second vector whose L elements are 
stored symbols of the transmitted training or estimated sequence ({*„}*{**„}). 
respectively, assuming a duration of the impulse response equal to L symbols; 
being Q(n) a matrix transformation of unitary rotation, implemented through a 

sequence of elementary transformations Q 0 ,Qj ,A , that applied to 

produces the cancellation of column g n introducing a zero for each elementary 

transformation; 

d5) determination of a weight vector h„ through the following expression: 
L s (n)V s (n)h n = h c t (n - l) , according to which said vector of optimum samples is 



opt ' 



15 updated as follows h (n) = , 

h n h n 

d6) updating of said optimum weight vector through the following expression: 
L'(n)w o >) = U>)h op ,(n). 

7. Equalization method according to claim 6, characterized in that said 
expression calculated at step d5) if divided into two calculation phases, introducing an 

20 unknown intermediate vector u rt defined by the following first equation: 

V s (n)h n = u n , being u„ a solution of the following second equation: 
L f (rt)u„ = h opt (n-l) that replaced in its turn in the starting equation gives 

simplifying the relevant calculation. 

8. Equalization method according to claim 6 or 7, characterized in that in said 

25 step e) said matrix L(« - l) is multiplied by VX, where X is an arbitrary real constant, 
lower than one unit, selected in relation to the non stationariness of the channel making 
said adaptive equalizationmore effective. 
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9. Equalization method according to any of the previous claims, characterized 
in that said estimate of the symbols of the transmission burst made in step d) is an 
estimate of sequence ({3„}) at maximum likelihood made through the Viterbi 
algorithm, using for this said optimized samples of the channel impulse response 

5 (h op ,(w)) to calculate the values of appropriate branch metrics of a trellis indicating all 
the possible status transitions of the receiver, the expression of the branch metrics 
being known in said receivers. 

10. Equalization method according to claims 1 and 9, characterized in that in 

the updating step e) the symbols of the training sequence ({*„}) are replaced by a low 

10 delay estimated sequence of symbols [{&„}) forcing the Viterbi algorithm (STIM) to 
emit decisions at the end of a time lower than that taken to decide said maximum 
likelihood estimated sequence [{B n }) . 




11. Equalization method according to claim 10, characterized in. that, 
coinciding with any new symbol of said low delay estimated sequence theL 

15 elements of said vector x{n) are replaced by a same number of symbols copied as a 
whole from the path selected in the trellis at the present symbol time, preventing by this 
the propagation of errors in the joint estimate of said optimum weights and channel 

impulse response samples (w 0/7/ (V)> h op ,(*))- 



12. Equalization method according to claim 10 or 1 1 , characterized in that said 



20 vector y(n) suffers a delay equal to that suffered as a whole by said low delay 



13. Equalization method according to any of the previous claims, characterized 
in that it includes also a step for the continuous monitoring of the power received on 
the antennas (ARY) and detection of possible peaks exceeding a given threshold. 

14. Adaptive equalization method according to claims 6 and 13, characterized 
in that, at an instant n of detection of said peak of power received on the antennas the 
following actions are accomplished: 

- advanced updating, within the symbol time n , of the sole triangular matrix L of said 
matrix L(w) using the vectors y(«) obtained from the samples of signals on the 
antennas, present in the stored bursts, starting from instant n up to an instant 




estimated sequence ( [3 n 
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n + AT that limits a transient for the setting of values calculated in the updating of 
said matrix L(«) ; 

- updating, for n < n < n + AT , of matrices U(«) and 'L,(/7) of said triangular 
matrix L(«) with said updating vector g„ and execution of said step d6) which 

5 updates said vector of optimum weights (w 0/w ), while keeping constant said 

samples of the channel impulse response (h op/ ) with values assumed at said 
instant »: 

- restoration, when time n + AT elapses, of the updating modes of said matrix L(n) 
preceding said instant n . 

10 15 Equalization method according to claim 14, characterized in that said 

advanced updating of the sole triangular matrix L , is made applying the generic said 
elementary transformation Q, to the first M elements of the Mh column of said 

triangular matrix L(«), where Mis the order of matrix L, repeating the elementary 

transformation Q, for all the values of index i and storing the transformations 

15 occurred. 

16. Equalizer in receivers of base stations of a cellular mobile system, in 
particular for the cancellation of isofrequential interferes in signals received by array 
antennas (ARY) of base stations of a cellular mobile system, said signals (y(«)) being 
received as symbol bursts of a transmission burst, including a training sequence 
20 ({**}) known to the receiver, said bursts received being demodulated, digitized, 
stored, and synchronized through discrimination of the training sequence, 
characterized in that it includes: 

a) a spatial filter (BMF) of signals coming from the antennas (y(«)) that combines 
among them corresponding symbols of the different bursts stored on the basis of 

25 appropriate weights (w op ,-(n)) ; 

b) a time matched filter (FADAT) , filtering the signal (y out {n)) coming from said spatial 
filter (BMF); 

c) means for the joint optimization (SYSTOL1 ,3,4,5) of said weights and samples. 
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(w op ,(w), h 0/ , r («)) on the basis of the symbols of said training sequence ({*„}) and 
of said signals coming from the antennas (y(n)) present in the stored bursts; 

d) decider (STIM) that estimates the symbols ({4}) of the transmission burst 

corresponding to the sequence coming out from the time matched filter; 

e) updating means (SYSTOL2, 3,4,5) of optimized values of said weights and samples 

( w o P /(' ? )> h o^(")) acting out of the time slot of the training sequence ({*„}) and 
including the previous means a)...d) in which the optimization means c) act now on 
the basis of symbols estimated ({£„}) by the decider d), making by this adaptive 
said equalization. 

17. Equalizer according to claim 16, characterized in that said joint 
optimization means (SYSTOL1 ,3,4,5) include a first systolic array (SYSTOL1) 
receiving said symbols of the stored bursts, organized as in an algebraic vector y(w) of 
M elements, as many as are the antennas, and the symbols of said training sequence 
{{ x n}) organized as in an algebraic vector x(w) of L elements, equal to the length of 

15 the channel impulse response, and generating, coinciding with each symbol time «-th, 
a matrix updating of the triangular matrix L(«) whose elements allow the updating of 
the estimate of a vector h («) storing said optimum samples of the impulse response 
and the estimate of a vector w opf (#?) that stored said optimum weights. 

18. Equalizer according to claim 17, characterized in that said vectors y(») 
20 and x(n) reach in column, in the order indicated, the input of the cells of said first 

systolic array (SYSTOL1 ) corresponding to the elements of the first column of said 
matrix h(n) , the whole of cells of the triangular matrix including the calculation means 
of the following matrix relation: 

" 1 P 



[a out] <- [a in] 



25 where: 

- a pointed by the arrow indicates the new content of the cell; 

- out is an output signal; 

- in is a signal entering each cell; and 
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p is a both input and output signal, calculated as follows: 



in 

P = -; 



the cells corresponding to the elements of the main diagonal of said matrix L(n) 
containing the alternative calculation means, such for which the signal p is only an 
output one, the signal out is absent, and the content a is updated as follows: 



a <— ^J\a\ 2 + l //7 i 2 • 

19. Equalizer according to claims 16 and 18, characterized in that said 
updating means (SYSTOL2, 3,4,5) include a second systolic array (SYSTOL2) whose 
cell include calculation means differing from those of said first systolic array 

10 (SYSTOL1) only due to the fact that the content a is multiplied by a constant numeric 

value VX , where A is an arbitrary real constant, lower than one unit selected in 
relation to the non stationariness of the channel to render said adaptive equalization 
more effective. 

20. Equalizer according to claim 18 or 19, characterized in that said triangular 
15 matrix L(«), image of said first and second systolic arrays (SYSTOL1 , SYSTOL2), is 

defined as follows: 

[LM 0 

L(n) = Low L.W. 

in which the elements L(w) and L,(«) are also triangular matrices, and U(w) is 
rectangular, or square. 

20 21. Equalizer according to claim 20, characterized in that said optimization 

means (SYSTOL1 ,3,4,5), or said updating means (SYSTOL2,3,4,5), include also a 
third systolic array (AUTOVET) and a fourth one (SYSTOL5) placed in cascade, said 
third array receiving said element L s {n) and udating the elements of said vector 
h opt (n) , and said fourth systolic array receiving said elements L{n) and U(w) and the 

25 vector h opt {n) , and updating the elements of said vector w opt {n) . 

22. Equalizer according to any claim 16 through 21, characterized in that said 
detector (STIM) is a Viterbi estimator performing a maximum likelihood estimation of 

the sequence of symbols ({&„}) of the transmission burst, using to this purpose said 
optimized samples of the channel impulse response h opt (n) to calculate the values of 
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appropriate branch metrics of a trellis that indicated all possible status transitions of the 
receiver, the expression of the branch metrics being known in said receivers. 

23. Equalizer according to claims 16 and 22, characterized in that it also 
includes: 

5 - storage means (DECTEMP) of the symbols belonging to an estimate path selected 
in the trellis, with a delay having duration AT lower than the emission delay of the 
symbols of said maximum likelihood estimated sequence the symbols 

stored being organized like an algebraic vector k'{n) reaching the input of said 
second systolic array (SYSTOL2); 

10 - delay means (RIT) that delay by AT said algebraic vector y(/i) before it reaches 

said second systolic array (SYSTOL2); 

- logic means (MID) for detection of coincidence, or not, between the position in the 
burst of symbols received and the time slot of the training sequence in symbols 

transmitted ({•*„}); 

15 - selection means (SEL), controlled by said logic detection means (MID), which, in 
case of coincidence with the training sequence, extend the bus where said group of 
values L(«) runs at the input of said third systolic array (AUTOVET) and, in parallel, 
at the input of said fourth systolic array (SYSTOL5). while in case of non- 
coincidence extend on the contrary the buffer transporting said group of values 

20 L'(n); 

- transcription means (CONT) of the contenut of said first systolic array (SYSTOL1 ) in 
the corresponding cells of said second systolic array (SYSTOL2), controlled by said 
logic means (MID) for the detection of the training sequence, that enable 
transcription simultaneously to the closing transition of the time slot of the training 

25 sequence. 

24. Equalizer according to claim 23, characterized in that it also includes: 

- detection means of an instant n of the occurrence of a power peak in said signals 
received on the antennas, and generating a signal of presence of strong interferer 
(PIC), active out of the time slot of the training sequence for said duration AT ; 

30 - stop means (SWT) of said extension towards the input of said third systolic array 
(AUTOVET), controlled by said signal of presence of strong interferer (PIC); 

- freezing means of the content (h op ,(n)) of said fourth systolic array (SYSTOL5) at 
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said instant n , controlled by said signal of presence of strong interferer (PIC); 
- means (COM, FIFO) for the variation of the internal configuration of said second 
systolic array (SYSTOL2), controlled by said signal of presence of strong interferer 
(PIC). 

5 25. Equalizer according to claim 24, characterized in that said variation means 

(COM, FIFO) of the internal configuration of said second systolic array (SYSTOL2) 
include electronic switches (COM) switching the lines of the array connecting the cells 
corresponding to the elements of the last line of said triangular matrix L(«) with the 
elements of the first line of said rectangular matrix U(«) , and connecting a same 

10 number of FIFO memories between the two ends of the switched lines. 

26. Equalizer according to any claim 17 through 25, characterized in that said 
first (SYSTOL1) and second (SYSTOL2) systolic arrays are grouped in a unique array 
having the common triangular structure of cells of said first and second array, and that 
it includes the additional multiplication means of said content a of the single cells for 

15 said constant numeric value VX, said multiplication means being controlled by said 
logic detection means (MID) of the time slot of the training sequence. 

27. Equalizer according to claim 21, characterized in that said third systolic 
array (AUTOVET), is divided in the cascade of two systolic arrays (SYSTOL3, 
SYSTOL4) having identical triangular structure, of which the upstream array 

20 (SYSTOL3) receives said group of values L,(«) , and the previous vector h opt {n- l) , 
of the downstream array, and calculates an intermediate vector u„ supplying it to the 
downstream array (SYSTOL4), reached also by a group of values L",(n), for the 
calculation of said vector h op ,( w ) • 

28. Equalizer according to any claim 16 through 27, characterized in that when 
25 said training sequence belongs to a short midamble versus the length of the whole 

burst, said equalizer is duplicated, the two equalizers processing in parallel one or the 
other semiburst at the sides of the midamble, and giving a relevant estimated 
subsequence to a burst reconstruction circuit which places in the right time order the 
symbols of the two estimated subsequences, obtaining a unique maximum likelihood 

30 sequence estimate of the data transmitted {d n } . 
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